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Abstract 

Kinetically constrained lattice models of glasses introduced by Kob 
and Andersen (KA) are analyzed. It is proved that only two behaviors 
are possible on hypercubic lattices: either ergodicity at all densities or 
trivial non-ergodicity, depending on the constraint parameter and the di- 
mensionality. But in the ergodic cases, the dynamics is shown to be in- 
trinsically cooperative at high densities giving rise to glassy dynamics as 
observed in simulations. The cooperativity is characterized by two length 
scales whose behavior controls finite-size effects: these are essential for 
interpreting simulations. In contrast to hypercubic lattices, on Bethe lat- 
tices KA models undergo a dynamical (jamming) phase transition at a 
critical density: this is characterized by diverging time and length scales 
and a discontinuous jump in the long-time limit of the density autocorre- 
lation function. By analyzing generalized Bethe lattices (with loops) that 
interpolate between hypercubic lattices and standard Bethe lattices, the 
crossover between the dynamical transition that exists on these lattices 
and its absence in the hypercubic lattice limit is explored. Contact with 
earlier results are made via analysis of the related Fredrickson- Andersen 
models, followed by brief discussions of universality, of other approaches 
to glass transitions, and of some issues relevant for experiments. 

1 Introduction 

The glass transition — whether a true transition or a sharp crossover in the 
dynamics — is one of the longest standing unsolved puzzles in condensed mat- 
ter physics. Many liquids, when cooled fast enough to avoid crystallization, 
initially remain liquid but at lower temperatures appear to freeze into solid-like 
structures that nevertheless show no signs of crystalline order. The time scales 
for structural relaxation in such metastable super- cooled regimes increase dra- 
matically as the temperature is lowered and materials are conventionally called 
glasses when the structural relaxation times are longer than those accessible 
in typical experiments. In fragile liquids the relaxation time r increases more 
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rapidly than Arrhenius [1,2] and can be characterized by an effective activation 
free energy d\ogr / d(l/T) that increases by as much as a factor of twenty over 
a narrow temperature range. [2] Moreover, near glass "transitions", relaxation 
processes typically become complicated [1], often characterized by stretched ex- 
ponential decay of temporal correlations. Concomitantly, persistent spatial het- 
erogeneities appear in both experiments [3,4] and numerical simulations [5-8]. 

Yet, in spite of a great deal of theoretical effort over the past few decades, 
a real understanding of these and related phenomena is still lacking. Indeed, 
the most basic issues are still unresolved: Is the rapid slowing down due to 
proximity to an equilibrium phase transition of some kind (albeit perhaps in a 
restricted part of phase space)? Or is the underlying cause entirely dynamical? 
The latter scenario is supported by the lack of evidence for a diverging static 
correlation length and the weak temperature dependence of structural correla- 
tions, [2] while growing dynamical correlation lengths have been found in both 
experiments and simulations. [3-8] Whether or not there is an actual transition 
of some kind, what is the physical mechanism responsible for the slowing down 
of the dynamics and the growth of dynamical length scales? 

Theoretical progress has been hampered by a shortage of models that capture 
at least some of the features believed to be essential, yet are simple enough 
to analyze. In other fields of statistical physics simple models have played 
vital roles leading to new understanding, to sharpening of questions, and to 
development of improved models. 

1.1 Kinetically constrained lattice models 

One of the few types of simple models for glasses are the kinetically constrained 
models, in particular, lattice models with constrained dynamics; for reviews see 
[9,10] . In this paper, we focus on a particularly simple class of models introduced 
by Kob and Andersen (KA) [11]. Before introducing these KA models and 
presenting our results, we give a short overview of kinetically constrained models 
more generally, including their motivation and the key questions about their 
behavior that could provide a deeper understanding of glass transitions in real 
systems. 

The basic Ansatz that motivates kinetically constrained models (see e.g. [12]) 
is that the key to glass transitions is the geometrical constraints on rearrange- 
ments of the atoms or molecules and the effects of these constraints on the dy- 
namics; static correlations, beyond those present in dense liquids, are assumed 
to play no role. The simplest models, which include KA models, are stochastic 
lattice gases with no static interactions beyond a hard core which restricts each 
site to be occupied by at most a single particle. The non-trivial aspect of these 
models arises from local constraints that are imposed on the motion of particles. 
Typically, the dynamics is given by a continuous time Markov process with each 
particle attempting, at a fixed rate, to jump to a randomly chosen empty neigh- 
boring site, but such a move is allowed only if the local configuration satisfies 
one or more constraints that depend on occupations of other nearby sites. 

The motivation for the kinetic constraints in conservative models arises from 
the behavior of molecules in dense liquids: the presence of surrounding particles 
can severely inhibit the motion of a molecule. In particular, a molecule can be 
caged by its neighbors in such a way as to stop it moving a substantial distance 
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until the cage is opened by the motion of other particles. When caging is ubiq- 
uitous, such local constraints can produce a degree of cooperative behavior that 
slows down the dynamics: one can readily imagine this to be the underlying 
mechanism that causes glass transitions. Since the idea behind the caging pic- 
ture is that glass transitions are purely dynamical phenomena, with any changes 
in static correlations playing a minor role (see e.g. [12]), it is natural to explore 
models whose dynamics satisfy detailed balance with respect to some simple en- 
semble, in particular the trivial measure that is uniform over all configurations 
with the same number of particles: this corresponds to ignoring static interac- 
tions beyond hard cores. [These models are therefore quite different from the 
statically constrained "lattice glass models" introduced in [13]]. Because of the 
neglect of interactions beyond hard cores, no equilibrium transition can take 
place in purely kinetically constrained models; clearly they cannot approximate 
the behavior of glass- forming liquids in all regimes, including the supercooling 
that avoids crystallization. Nevertheless, despite their simplicity and discrete 
character, kinetically constrained lattice models might still capture, at least on 
some range of time scales, the key dynamical aspects of real glass transitions. 
Indeed, numerical simulations show that some such models do display sluggish 
dynamics and much phenomenology that is reminiscent of glassy behavior and 
its onset as density increases [9, 10]. There has been a renewal of interest in 
these models for explaining the behavior of glass-forming liquids [12,14-24,40]. 
Kinetically constrained lattice gases have also been studied in the context of 
granular systems [25] which also display a glass-like dynamical arrest as their 
density is increased; this so-called jamming transition [26,27] occurs at a density 
well below close-packing density and can thus not have purely entropic origin. 

Note that there is another class of kinetically constrained models whose ele- 
mentary degrees of freedom are spins and not particles. An important example 
on which we will focus on sectionals provided by the Frcdrickson and Andersen 
(FA) models [28]. These are the non-conservative version of the previous ones: 
there is no static interaction among spins on different sites and the non-trivial 
aspect arises from the local constraints imposed on the possible spin flips. In this 
case, the motivation for kinetic constraints comes from the concept of dynamic 
facilitation [12,28]: the evidence of dynamical heterogeneities, i.e. the existence 
regions of space with very different relaxation times suggests that, on a coarse 
grained level, supercooled liquids might be considered as a mixture of mobile 
and non-mobile regions. The non-conservative character comes from the fact 
that the ratio of mobile to non-mobile regions is non-conserved in time. Kinetic 
constraints arise from the fact that the vicinity of mobile regions is required in 
order to enable (facilitate) mobility in a non- mobile regions. 

1.2 Characterization of possible glass transitions 

The most crucial questions about kinetically constrained models of glasses are 
whether dynamical arrest takes place with increasing density, and, if so, whether 
this is the result of some type of sharp transition or instead a gradual freezing. 
In either case, one needs to understand the cooperative nature of the dynam- 
ics: in particular, the emergence of growing length scales that characterize the 
coopcrativity. 

We briefly sketch several possible scenarios for glass transitions. The most 
striking possibility would be an actual ergodic / non-ergodic transition at some 
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critical density: this possibility, which corresponds to an analogous conjecture 
for the nature of the dynamical arrest in glass forming liquids, is the primary 
subject of the present paper. Because of the simple nature of the KA models, 
this question can be completely separated from possible changes in equilibrium 
static behavior. As mentioned above, for kinetically constrained models whose 
dynamics conserves the number of particles and satisfies detailed balance with 
respect to the uniform (non-interacting) measure, the trivial distribution that 
is flat over all configurations with the same density is always stationary and 
there can be no equilibrium transition. Nevertheless, due to the constraints on 
the allowed particle moves, ergodicity could be broken. Non-crgodicity means 
that even in the limit of long times, time averages of local quantities would 
not converge to averages over the static equilibrium ensemble. This certainly 
occurs for some kinetically constrained models on finite lattices for which the 
configuration space typically breaks into disconnected irreducible components: 
i.e. there exist two or more sets of configurations that can not be connected 
to one another by any sequence of allowed moves. A key question is whether 
this non-ergodicity can persist in the thermodynamic limit for which the grand 
canonical ensemble at fugacity p (Bernoulli product measure at density p) is 
an equilibrium distribution. Specifically, can an ergodic/non-ergodic transition 
occur at some density p c such that the dynamics on the infinite lattice is ergodic 
for p < p c , but not ergodic for p > p c ? Such an ergodicity-breaking transition 
is often considered an ideal glass transition; it has been advocated as the pri- 
mary cause of the slowing down of glass-forming liquids. Ergodicity breaking 
transitions have been found to occur in several approximate theories, includ- 
ing in mode coupling theory [29] and in fully-connected quenched random spin 
models [67]. 

A weaker type of dynamical arrest could in principle take place even if er- 
godicity is not broken: in particular, a diffusive/ sub- diffusive transition of the 
dynamics. Indeed, the motion of a tracer particle could become non-diffusive at 
high densities: this could be observed by singling out of an infinite system that 
is in equilibrium, one particle, the tracer, and following its motion. There exists 
a general rigorous result that restricts the behavior: as long as the dynamics is 
ergodic, the position of a tagged particle converges under diffusive space-time 
rescaling to Brownian motion with a density-dependent self-diffusion coefficient 
Ds{p) > [30,31]. But in order for the tagged particle to diffuse, the self- 
diffusion coefficient has to be non-zero: for the normal hard core lattice gas 
(simple symmetric exclusion), this has been proven [30] and in [32] the result 
has been extended to all the models such that the rate of any given jump to an 
empty neighbor is bounded away from zero. For kinetically constrained lattice 
gases, however, the tagged particle could be slowed enough that D$ becomes 
zero at sufficiently high density. A diffusivc/sub-diffusive transition would then 
presumably take place at some critical density, po with D$ > for p < pu but 
Ds = for p > po- Such a transition has been conjectured to occur in real 
glasses: and several experiments on supercooled liquids near their glass transi- 
tion indeed show a striking decrease of the mean square displacement of a tracer 
particle over the observed range of times [2]. 

Another way in which diffusion could potentially break down involves the 
breakdown of the conventional hydrodynamical limit that usually holds on long 
length and time scales [31]. For lattice gases (and more generally for systems 
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with particle number conservation but no momentum conservation) the coarse- 
grained density profile evolves, at long times, via the diffusion equation, dtp = 
V(D(p)Vp), with the diffusivity, D(p), dependent on the local density. For 
generic kinetically constrained stochastic models, proving the existence of the 
hydrodynamic limit, as has been done for conventional hard core models, is not 
a trivial task. The natural conjecture, if a transition of some kind did occur, 
would be the vanishing of the macroscopic diffusivity, D(p), at high densities. 
This would give rise to a form of macroscopic arrest, in which density profiles at 
high densities would evolve, if at all, only sub-diffusively. Note that, in contrast 
to the diffusion of a tagged particle, evolution of density profiles need not involve 
motion of individual particles far from their initial neighbors: density diffusion 
and self diffusion thus probe rather different aspects of the dynamics. 

Consideration of the types of macroscopic diffusion raises the more general 
question of relaxation towards equilibrium: this could be anomalously slow even 
if the system remains ergodic at all densities. An interesting quantity to consider 
is the dynamical structure factor, the Fourier transform of the density-density 
correlation function at non-zero wave vector. For normal lattice gases, the dy- 
namical structure factor decays exponentially at long times but for kinetically 
constrained models, the relaxation could perhaps be slower than exponential, 
at least above some critical density. Indeed in numerical simulations of some 
kinetically constrained models at high density, sub-exponential relaxation has 
been found over a range of more than three decades [11]. This, as well as 
the relaxation of different response functions (see e.g. [14,33]) is usually fitted 
with a stretched exponential S(t) ~ exp[— (t/r)^], with /3 < 1 (also called a 
Kohlrausch- Williams- Watts function). Many glass- forming liquids show similar 
behavior with relaxation functions well-fitted by simple exponentials at high 
temperatures, but by stretched exponentials at low temperatures [2,34,35]. In 
glass-forming liquids this behavior arises from the superposition of many relax- 
ations with different decay rates, each corresponding to different spatial regions 
within a dynamically heterogeneous spatial structure which persists for longer 
time scales. [3]. If such sub-exponential relaxation does occur in kinetically 
constrained models — even if only over a limited range of times — a detailed 

understanding of it would shed light on the behavior of more realistic systems, 
l 

Numerical simulations of Kob-Andersen models, as reviewed below, have 
found very sluggish dynamics at high density [11], with much phenomenology 
that is very similar to glass-forming liquids. Is this a signature of some type of 
real dynamical glass transition at a critical density? More generally, does one 
exists for at least some versions of the KA models? If not, what is the basic 
mechanism for the extremely rapid onset of the slowing down as the density is 
increased? Are there collective processes involved? What if any, are the char- 
acteristic length scales, and how are the time scales related to these? What 
is the physical mechanism giving rise to the dynamical heterogeneities? And, 
most importantly, what can one learn from the KA models about the behavior 
of more realistic models? In this paper we will address all of these, although 
the last, regrettably, only to a limited extent. 



It is important to note that stretched exponential relaxation can occur in the absence of 
such complications as heterogeneities: indeed, the two dimensional Ising model with non-spin- 
conserving dynamics exhibits just such behavior in its ordered phase. [36] 
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Some of the results of this paper have already been presented in a short 
letter [37] and the proof of the positivity of the self-diffusion coefficient for KA 
models was presented in [38] ; these both made use of the results derived in the 
present paper. 

1.3 KA models: definition and heuristics 

Kob- Andersen models were introduced [11], as mentioned above, to test the 
conjecture that cage effects in liquids can induce a sudden onset of dynamical 
arrest as the density increases thereby being responsible for a glass transition 
in such systems. The conditions for particles to be able to move in KA mod- 
els mimic those caused by the geometrical constraints imposed by surrounding 
particles on the motion of molecules in dense liquids. 

The Kob- Andersen models we consider are natural generalizations of those 
of reference [11]: they are defined on d— dimensional hypercubic lattices A G Z d 
with zero or one particles allowed per lattice site and a parameter, to, that 
represents the set of allowed particle motions as explained below (for a more 
formal definition see section T2.5J1 . Each particle attempts, at rate unity, to move 
to a randomly chosen one of its z = 2d nearest neighbor sites. But a move 
from a site i to a neighboring site y is allowed only if site y is empty and 
the particle has no more than to occupied neighbors both before and after the 
move. If any of these constraints is not satisfied the particle remains at site x 
(at least until its next attempt to move). The possible values of to range from 
m = 0, . . . , 2d — 1, the upper limit being one less than the coordination number, 
z (due to the hard core constraint). In finite systems, one can consider a fixed 
number, N , of particles. But in the infinite systems of primary interest, the 
particle number density p, must be specified instead. 

As the interesting behavior occurs at high densities — p close to one — it is 
useful to reformulate the rules in terms of vacancy motion, i.c dynamics of the 
empty sites. As can readily be verified, the above constraints with parameter to 
correspond to vacancies allowed to move only if both the initial site before the 
move, and the final site after the move, have at least s = z — to — 1 neighboring 
vacancies. The hypercubic KA models are thus completely defined by the pa- 
rameters d and to, or equivalently, d and s; we will usually use the latter. The 
special case m = 2d — 1 corresponds to s = 0: this KA model thus corresponds 
exactly to conventional hard-core lattice gas dynamics (i.e. simple symmetric 
exclusion process). 

For all values of s the equilibrium static properties are identical to those 
of the hard-core lattice gas model with no other interactions, and are hence 
trivial. But the dynamics depends radically on s. For s > d, the dynamics is 
so constrained that KA models are trivially non-ergodic at any p £ (0, 1]. This 
can be checked by noticing that any fully occupied hypercube of can never be 
broken up, therefore at any finite density there exists a finite fraction of particles 
which are forever blocked. We thus focus on < s < d. 

Why might one expect interesting dynamics for s > 1? In contrast to the 
conventional lattice gas (s = 0) for which individual vacancies can move freely, 
for non-zero s on hypercubic lattices, there are no finite mobile clusters of va- 
cancies that can move freely — albeit together — in an otherwise fully occupied 
system. By this we mean that, if the surrounding lattice is completely filled, it 
is not possible to construct a a path of allowed nearest neighbor moves through 
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which one can shift the cluster of vacancies. This follows from the existence of 
infinite frozen configurations which, as none of the particles in them can move at 
all, impede the motion of any cluster of vacancies. The simplest infinite frozen 
configuration is a fully occupied slab of particles that is infinite in d — s of the 
lattice directions and has width two in the other s directions: by construction, 
none of the particles in such a slab can move with the KA dynamical rules 
with parameter s. [In contrast to the absence of mobile clusters of vacancies 
on hypercubic lattices with non-zero s, on a two-dimensional triangular lattice 
with s = 1, isolated vacancies cannot move but neighboring pairs of vacancies 
can move freely and equilibrate the system (see [38] for a detailed and rigorous 
discussion).] 

Although for non-zero s on hypercubic lattices there are no finite mobile clus- 
ters of vacancies, neither are there finite frozen clusters for s < d. This follows 
because a particle at a corner of any fully occupied region that is surrounded 
by empty sites can move; thus on an infinite lattice there are no finite frozen 
clusters. Nevertheless, there do exist infinite frozen sets of particles such as the 
completely filled slab defined above. Such simple frozen configurations have very 
little entropy — in fact, zero entropy density — but a much larger set of frozen 
configurations exist that does have non-zero entropy density. For example, for 
the simplest case, s = 1 on a square lattice, in addition to configurations with a 
single two-wide frozen slab that is infinitely long, any configuration with a set 
of such slabs that each begin and end at T-junctions with other such slabs will 
be frozen. A crucial question about static configurations is thus whether, at 
sufficiently high densities, almost all infinite configurations contain an infinite 
- and hence connected — set of frozen particles. As we shall see, answering 
this is by no means trivial. For all the nontrivial cases, the above observations 
imply that the sluggish dynamics at high densities must be intrinsically collec- 
tive, involving rearrangements of increasingly larger and larger regions as the 
density increases. 

1.4 Previous numerical results 

Before introducing the analytic methods we will use, we briefly summarize the 
results of various numerical simulations of the best studied KA model: the 
three-dimensional case with s = 2. It has been suggested that these simulations 
provide support for the conjecture of an crgodic / non-ergodic transition at a 
non-trivial density, p. In reference [11] the self-diffusion coefficient of a tagged- 
particle was measured and the results found to fit well a power law form that 
vanished at a finite density: D s cx (p — p) a with a ~ 3.1 and p ~ 0.881. 
Approaching this same apparent critical density the rate of temporal relaxation 
of density-density correlations appeared to vanish as an inverse power of p — p. 
Both results — which hold over roughly four decades of Ds and two decades of 
p — p — were strongly suggestive of a dynamical glass transition at p = p. 

Furthermore, the asymptotic decay of the Fourier transform of density- 
density correlation function was fitted with the form S(t) ~ exp[(— £/ r ) ] with 
an exponent j3 close to one (exponential decay) for low and intermediate den- 
sities and decreasing monotonically with p (stretched exponential decay) at 
sufficiently high density (for p > 0.75). 

Later work found that other typical features of "glassy dynamics" occur at 
high densities in this d = 3, s = 2 model. 
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In [20] numerical simulations found dynamical heterogeneities to be persis- 
tent for very long times, at densities higher than roughly p. In [21] the effects of 
boundary sources of particle — corresponding to a grand canonical ensemble — 
were studied, in particular by quenching the chemical potential p of the sources 
below the apparent critical value that corresponds to the equilibrium fx(p). A 
non-equilibrium regime was evident, and "aging" effects similar to those ob- 
served in glasses were found. Properties that play roles in various theoretical 
approaches to glass transitions have also been studied. In particular, the rela- 
tionship between the "configurational entropy" of the number of distinct frozen 
configurations and an "effective temperature" which emerges from the aging 
dynamics was investigated to test Edward's hypothesis [22] as formulated by 
Kurchan in the context of aging dynamics; this had been found in certain mean 
field models studied in the context of glasses, see [23]. Edwards' hypothesis 
appears consistent with numerical results on the d = 3, s = 2 KA model [22]. 

The results announced in [37] and explained in detail in the following provide 
a theoretical framework to interpret numerical simulations. Two recent works 
[14, 24] have indeed analyzed the d = 3, s = 2 KA model taking into account 
our results and checking some of our analytic predictions. 

1.5 Summary of exact results 

From the numerical simulations described above, it is clear that the dynamics of 
some KA models is very sluggish at high density and displays much of the basic 
phenomenology of glassy dynamics observed in many experimental systems. 
Furthermore, there appears to be a rather sharp onset of the glassy behavior 
as the density is increased. The analytical methods of this paper yield an 
understanding of the mechanisms underlying the slow dynamics of KA models 
and address the basic qualitative and quantitative questions about them. 

The primary question is whether there exists some kind of actual dynamical 
transition for the d — 3 s = 2 (or other) KA model, or whether the apparent 
"transition" seen in simulations is only a result of the finite size and finite time 
limitations of such simulations. We first investigate whether an ergodic/non- 
ergodic transition can occur at a non-trivial density in KA models. We conclude 
that, for any dimension and any choice of the constraint parameter s < d, 
an ergodic / non-ergodic transition cannot take place at a non-trivial density 
p € (0, 1). In other words, the sluggish behaviour found in simulations cannot be 
the mark of an "ideal" glass transition. More quantitatively, we derive a density 
dependent characteristic length scale which separates two different regimes: 
for samples with linear dimensions L <C the configuration space breaks 

into many ergodic components; while for those with L 3> S(p) a single ergodic 
component dominates — although other disconnected components still exist (a 
precise definition of is given below). The length, S, whose dependence 
on p varies with d and s, diverges at high density: ~ expexp (c/(l — p)) 

for the d = 3 and s = 2 model. Therefore, large systems of size L 3 , are likely 
to appear ergodic only if p < p c (L) = 1 — (c/ log log L). Although for smaller 
sizes the apparent ergodicity breaking takes place at lower values of the density, 
the dependence on L is very weak for large L. Of course, to properly analyze 
data from simulations, knowledge of the existence of a cross over length is a key 
ingredient for disentangling the finite size effects which will always be limiting 
at high densities. 
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The fact that ergodicity holds in infinite systems at all densities for s < d, 
does not rule out the possibility of a more subtle type of dynamic transition 
at a non-trivial density. In particular, a diffusive/sub-diffusive transition, if it 
existed, would explain the apparent vanishing of the self diffusion coefficient 
found in the above-mentioned numerical simulations [11]. We have studied the 
asymptotic behaviour of the diffusion coefficient, D$, of a tagged particle. The 
result, discussed in detail in [37] and [38], is that D$ is strictly positive at 
any density p < 1, so that a diffusive / sub-diffusive transition cannot take 
place. Moreover, this analysis [37] unveils the nature of the collective processes 
which are needed for both equilibration and diffusion at high density. The 
characteristic time scale, r, of the slow cooperative dynamics can be calculated 
and for the non-trivial cases, namely s < d, r diverges faster than any inverse 
power of 1 — p as p — ► 1. 

Given these results for general KA models on hypercubic lattices, an im- 
mediate question arises as to whether KA models can ever exhibit a transition 
at a non-trivial density: if this can occur, it raises the possibility that the 
surprisingly sharp onset of the slow dynamics observed in simulations in, for 
example, s = 2, d = 3, might be due to the "ghost" of such a transition. Expe- 
rience suggests that the most likely models in which to find actual transitions 
arc "mean- field-like" models [67]; thus we study KA models on tree structures, 
more precisely, Bethe lattices, which often provide realizations of mean field 
approximations. We find that on Bethe lattices there exists a critical density, 
p c , above which ergodicity is broken and phase space breaks up into many dis- 
connected components. The transition at p c has aspects of both first order and 
second order transitions. 

An interesting question that we leave open for future investigation is the 
hydrodynamic behavior of KA models on hypercubic lattices. Although we 
conjecture that on large enough length and time scales hydrodynamics will be 
valid at any density less than one for ergodic KA models, this might be very 
hard to prove. Assuming that the conjecture is correct, an important issue, 
related to both experiments [39] and simulations [8] of glass-forming liquids, is 
the characteristic length and time scale beyond which hydrodynamic behavior 
sets in, in particular, how these scales increase with density. 

A recent interesting analysis of this behavior within Kinetically Constrained 
models have been presented in [40] . 

1.6 Outline 

The organization of the remainder of this paper is as follows: In section El we 
analyze s < d KA models on hypercubic lattices and prove ergodicity for any 
density p < 1. First, in l2.1l wc analyze the simplest interesting case, d = 2, s = 1, 
and prove that with unit probability there exists an irreducible component of the 
configuration space in the thermodynamic limit and find an upper bound, S"(p), 
on the crossover length S(p). In 12.21 we extend these results to the case d = 3, 
s = 2 originally considered by Kob and Andersen, and in 12.31 derive analogous 
results for s = 1 in general d. The convergence of large systems to the infinite 
system behavior is studied in 12. 1.41 we prove that for L > the probability 

of the maximal irreducible component goes to one at least exponentially rapidly. 
Using these results, in 12.41 we extend the irreducibility proof and find upper 
bounds for the crossover length for all s < d. Finally, in 12.51 wc prove that 
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irreducibility of the configuration space implies ergodicity. 

Section|3]is devoted to quantitative results on the cross-over lengths. In l3.ll 
by using recent bootstrap percolation results to obtain lower bounds on S(p) and 
combining these with the upper bounds, we find the asymptotic behavior of the 
density-dependence of the cross over length, up to an undetermined constant. 
In section EQ1 we calculate the exact value of this constant for the case d = 2 
s = 1 via the identification of the dominant mechanism which restores ergodicity 
in large systems at high density. In 13.31 we formulate an analogous conjecture 
for the dominant high density behavior of the original KA model, d = 3 s = 2. 

In section 0] the physical picture for the cooperative high density dynamics 
(which was presented in [37]) is explained and related to the results derived in 
this paper; details are left for a forthcoming paper. 

In section[S]KA models on Bethe lattices are studied. After recalling (section 
15.1(1 the definition and some properties of Bethe lattices, in 15.21 and 15.31 recur- 
sive relations are derived whose fixed point solution implies that an ergodicity 
breaking transition takes place at a non-trivial critical density p c , for any s < k 
— the branching parameter of the Bethe lattice. The transition is shown to be 
discontinuous for s < k — 1 fsection l5.3.1fl but continuous for s = k — 1 (section 
I5.3.2fl . The existence and location of this transition are related to a boostrap 
percolation transition in l5.3.3l In section l5~^l the focus is on quantitative results 
for Bethe lattice models: wc analyze in some detail the character of the transi- 
tion in the cases k = 3, s = 1 and k = 5, s = 3 and present results of numerical 
simulations for the k = 3, s = 1 model. The behaviour of the density-density 
correlation function and the corresponding susceptibility are studied, along with 
the configurational entropy, S c (p). By establishing a lower bound we prove that 
S c (p) jumps to a non-zero value at p c . We end this section, in 15.51 by ex- 
tending to KA model on decorated Bethe lattices: graphs with finite size loops 
that interpolate between simple Bethe lattice and finite-dimensional hypcrcubic 
lattices. 

In section[l)]we analyze the related Frcdrickson Andersen (FA) models, prov- 
ing, in !6.2l ergodicity at any finite temperature (completing an almost complete 
proof that had been given by Fredrickson and Andersen themselves). In 16. 31 we 
give quantitative predictions on length and time scales for the d = 2 two spin 
facilitated FA model. 

Finally, in [71 we draw conclusions and discuss some possible connections of 
our results to other theoretical treatments of glass transitions. In particular, 
the analogy of KA models on Bethe lattices to fully connected random p-spin 
models. We also comment briefly on possible extensions of KA models and 
connections to experiments. 

2 Rigorous results for KA models on hypercubic- 
lattices 

In this section we analyze KA models on hypercubic lattices with the parameter 
s in the non-trivial range 1 < s < d— 1. We prove that for any density p < 1, on 
an infinite lattice such KA models are ergodic. In contrast, as explained earlier, 
for larger values of s, s > d, KA models on hypercubic lattices are non-ergodic 
for all positive p because of the existence of finite frozen clusters. There are 
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thus only two possible behaviors for hypercubic lattices: either the KA process 
in the thermodynamic limit is ergodic at any density, i.e. p c = 1, or else it is 
never ergodic, i.e. p c = 0; therefore ergodic /non-ergodic transitions cannot take 
place in KA models on hypercubic lattices for any d and s. 

The strategy of the proof is as follows. Let A 6 Z d be an hypercube of linear 
size L. First we identify a component of the configuration space fl\ = {0, 1}' A I 
and show that it is irreducible, i.e. any two different configurations belonging 
to this component can be connected to each other by a sequence of elementary 
moves allowed by the KA rules. Second, we prove that, with respect to the 
natural measure on the space of configurations, Bernoulli product measure HA,p, 
the probability of this single irreducible component goes to unity for L — > oo. 
Finally we prove that, thanks to the product form of Bernoulli measure, the 
existence of an irreducible component with unit probability implies ergodicity for 
the infinite system, A = Z d . 

In the first four subsections, we give the first and second steps of the proof 
for all choices of d and s, leaving to 12.51 the last step, which is proved by a 
general argument that is independent of the specific parameters. 

We start with the simplest non-trivial model, d — 2, s = 1, giving the most 
details for this case. 

2.1 Irreducibility for d=2, s=l 

The only non-trivial model on a square lattice is the case s = 1: vacancies can 
move if and only if they have a neighboring vacancy both before and after the 
move. 

2.1.1 Frameable configurations 

In order to construct the large irreducible component of configuration space 
we start by defining a subset of the configurations in it. Specifically, let us 
define framed configurations of any square or rectangle as those in which all 
the boundary sites are empty (see figure We then define as frameable any 
configuration from which, by an allowed sequence of elementary moves, a framed 
configuration can be reached. As we show below, any two framed configurations 
that have the same number of particles can be connected to one another by a 
sequence allowed moves; the same then applies to all frameable configurations 
by definition. 

2.1.2 Irreducibility of frameable set 

The irreducibility of a set of frameable configurations can be checked as follows. 
Note that even though it is not needed for the present purposes we shall prove 
irreducibility considering the particles as distinguishable. This result is impor- 
tant for the proof of the positivity of the self-diffusion coefficient [38] . Consider 
a pair of neighboring sites with for example particles A and B respectively 

on {i,j} and j immediately to the right of i. To prove the needed result, it is 
enough to show that, for any choice of the framed configuration, it is possible 
to perform the permutation of A with B leaving at the end of the process all 
the other particles in their initial configuration. 
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Figure 1: A 6 by 6 framed configuration. Filled dots are occupied sites; empty 
dots vacant sites. 

Starting from the bottom right corner it is possible to raise the bottom row 
of holes as shown in first line of figure|21 and likewise starting from the top right 
corner, to lower the top row of holes. This procedure can be iterated until the 
row that contains sites i and j is "sandwiched" between two rows of holes (first 
configuration on second line of figure [2J. At this point, using the surrounding 
vacancies, one can easily construct a path performing the desired permutation 
of particle A and B from site i to j (second and third line of figurcEJ. Following 
this exchange, the initial configuration of the rest of the lattice can be restored 
by reversing the moves of the two rows of holes after which all particles other 
than A and B have returned to their original sites. The full procedure can be 
performed if instead there is a particle and a hole on Since any framed 

configuration can be transformed to any other framed configuration with the 
same number of particles through a sequence of permutations of nearest neighbor 
particles and/or holes, all the framed configurations with the same number of 
particles belong to the same irreducible component. By definition, so also do 
all frameable configurations, a much larger set. 

2.1.3 Frameable probability 

The next step is proving that, of all configurations with density p, the fraction 
that are frameable, /za iP (J 7 ), approaches one in the limit L — > oo. We will do 
this by first constructing a large subset of J- '. 

Consider a four by four configuration which has at its center a two by two 
square of holes — a two by two framed configuration — and in the surrounding 
shell at least two holes adjacent to each side of the inner square. It is easy to 
check (see figure [3J that such a four by four configuration is frameable. This 
procedure can be iterated to grow an L by L frameable configuration starting 
from a two by two nucleus of vacancies and requiring at least two vacancies in 
each side of each subsequent shell. Therefore, ^a iP (J 7 ) is bounded from below 
by the probability, ii^ }P {T Q ), of frameable configurations constructed with this 
procedure from a two by two nucleus of vacancies centered at the origin: 

(L-2)/2 

MA,p(-n>MA, P (^°) = (i-p) 4 n (i-p'-v^i-p)) 4 . (i) 

1=1 

Note the as the size of the shell grows, the probability that there are the requisite 
number of vacancies on each side increases until becoming close to unity for 
sufficiently large shells: specifically when I ^> 1/(1 — p). The large L behaviour 
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Figure 2: Sequence of allowed elementary moves starting from a generic 6 by 
6 framed configuration, whose net result is the exchange of two neighbouring 
particles, A and B, inside the square. Figures in the first line show the ele- 
mentary moves connecting the initial configuration to a configuration with the 
bottom row of holes raised by one. By analogous moves (skipped in the figure) 
one can lower the top row of holes of two steps. Now (first figure of second row) 
particles A and B are "sandwiched" between two rows of holes. The remaining 
figures show how to exchange particles A and B starting from this configura- 
tion. Note that no other particles have been moved during this exchange. The 
configuration which is equal to the initial one (first line, first figure) on the sites 
not occupied by A and B can subsequently be restored by using sequences of 
moves analogous to those in the first line to raise the one full row of holes and 
lower the other one. 



of log (pa, p (^ )) is determined by 

log(l-p)+^log[l-p 2i -2«/- 1 (l- /J )] ~ log(l-p)-5> 2i + 2lp 2l -\l-p)] , 
l i 

(2) 

which is a converging series for any p < 1. Therefore p^ tP {T°) converges to a 
non-zero limit when L — > oo at fixed density. 

Because of the divergence of the above series for p = 1, additional care is 
required to analyze the behavior when both p — ► 1 and L — > oo. Let a(£) = 
log(l — p — 2£p 2e ~ 1 (l — p)). Since a(£) is increasing in £, the following inequality 
holds: 
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Figure 3: Growing 4 by 4 framcablc configurations from a 2 by 2 framed con- 
figuration 
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Figure 4: An 8 by 8 frameable configuration obtained from the growing proce- 
dure. The dashed square indicates the position of the two by two seed. Note 
that the density in such configurations is small near the seed but higher farther 
away. 



4a(l) + 4 / di a(£ < 
log/iA, P (^°)-41og(l-p) < 

4 / di a(£) (3) 



where the last term can be rewritten by using the change of variables y = p 21 
and expanding for p ~ 1 as 

/W)* ' r+ W-V + y*') (4) 



2(1 - p) Jo 

Therefore, by combining equations 10 and Q and using Um j0 _ > .i(l — p) log(l — 
p) = we get 

lim hm (1 - p) lo g , A , p (F°) = 2 f 1 ^ l0S(1 -^ + ^ 1 °^ ) (5) 
i.e., in the high density limit 2 

lim MA,p(^°) = Mco,p(^" ) =i e-(^) (6) 



L — >oo 

with 



Moreover from 



1 ^ log(l -y + y\ogy) 



2 Note that here and in the following we will use for simplicity the loose notation 
p,cc,p{3~°) — e although the correct mathematical expression is log (J.oo.piJ 70 ) — 



\ up to vanishing corrections in the p — » 1 



limit. 
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logMoo,p(^" ) - log pi A , pC? 7 ") ^ 



2 log(l-y + ylogy) 

i - p Ji y 

2 f pL 
/ dy(-l + logy) 

-H-(/log[/]-2p L ) 
1 — p 




(8) 



we find that piA.piJ 70 ) ~ Poa,p{F°) for L 3> £2,1 (/?)> denning a characteristic 
length, £ d ,s{p), with 

6,1 =-Iog(l-p)/(l-p) , (9) 

which will play an important role in what follows. 

We now need to get from the probability, p,A, P (^F°), of a frameable region 
centered at the origin, to the needed result: the probability that a square of 
size L is frameable about some center, p,A,p{^F), converges to unity for L — > 00. 
To do this one must consider all the possible positions for the initial nucleus 
used in the growing procedure 3 . Naively, we could simply make the argument 
that in a very large system, it is extremely unlikely that the system would not 
be frameable around at least one of the many possible nuclei. But additional 
work is required to turn this argument into a proof since the events that the 
whole square is frameable starting from different nuclei are not independent. 
We postpone the actual proof to the end of subsection 12.31 where we consider 
the generic d-dimensional case, and here just outline the argument. 

The key idea is the following: from the above definition of £, the conditions 
needed to further expand a frameable square of size L > ^ are satisfied with 
probability close to one. Therefore, whether or not one can grow frameable 
configurations starting from nuclei at the centers of distinct squares of size £ are 
almost independent. By considering that there are O (L 2 /£ 2 ) disjoint squares 
of size £ inside the L by L square, A, we conclude that p.A.pi^) is almost one for 
L ^S> £/ yj /J,oo,p(J~°) ~ exp [c/(l — p)] and approaches one as L — > 00. The rig- 
orous version of this argument given in 12 .31 completes the proof of irreducibility 
in the thermodynamic limit. 

Let us now turn to the analysis of finite size effects. As stated above, the 
probability of the irreducible frameable component goes to one when L — > 00 
at any fixed density p £ [0,1). On the other hand, for any finite and fixed 
value of -L, the probability of any irreducible component is strictly smaller than 
one: blocked configurations exist and therefore the configuration space is never 
covered by a single component. Furthermore, when p — ► 1 the probability 
pA,p(M) of the maximal irreducible component (the irreducible component with 
the highest probability) goes to zero since blocked configurations are more and 
more important and each one is in a different irreducible component. Therefore 
linip^i p\ tP (A4) = at fixed lattice size while limx^oo p\ :P (J\4) = 1 at fixed 
density. Thus the simultaneous L — > 00 and p — > 1 limit depends on the 

3 The results obtained hold for the KA model on a finite square lattice with periodic bound- 
ary conditions. In this geometry, any point can be taken as the origin of the growing procedure. 



15 



relationship between L and p as the limit is taken: this is typical of finite- 
size behavior near phase transitions. The crossover length S2.i(j°) IS defined 
to separate the "large" and "small" system behaviors so that by sending L — ► 
oo and p — ► 1 with L ^> H2,i(p) or L <C H2,i(p), respectively, one obtains 
p,\ iP (A4) — > 1 or /iA.p(A^) — ► 0, respectively. This defines the limiting behavior 
of H(/j) up to an overall constant factor. A more precise definition is possible, 
for example, by the sequence of {L,p) pairs for which PA,p(M.) — \- Above 
results on the probability of the frameable component establish therefore an 
upper bound on the crossover length: 



-2,1 (p) < -2,i (p) ~ exp 



(10) 



.(1-P). 

with c defined in Q. The reason why S" is only an upper bound for the 
cross over length is because while our framing argument guarantees that for 
L > S u a single irreducible component (the frameable one) dominates, it says 
nothing about the probability of the maximal irreducible component, which 
could already be close to unity for much smaller sizes. Indeed, if we have not 
got roughly the correct size-dependent conditions for frameability, Lx L squares 
could be almost ergodic even for (1 — p) logL <C (1 — p) log E". This possibility is 
ruled out in section lXTI bv establishing a lower bound S' that has similar density 
dependence to H", ensuring that we do indeed have the correct asymptotic form, 
albeit with the wrong coefficient c. The exact coefficient will be calculated in 
section EP1 thanks to the definition of a different framing procedure. 

The cross-over length, S, is also the scale above which frozen sets of particles 
are unlikely to exist. In particular, a spanning network of fully-occupied two- 
wide slabs that all start and end either at system boundaries or at T-junctions 
with other slabs — a frozen configuration — is likely to exist for squares of size 
L < 5, but unlikely to exist for much larger squares. How does this occur? 
In particular, if a number of squares of size somewhat smaller than 5 are put 
together to make one of size a few times 5, what happens to the (formerly) 
frozen bars that terminated inside one of the original squares? A crucial feature 
is the fragility of a typical frozen network of slabs: if one slab is cut, as will 
occur if it ends on an edge of one of the now-interior smaller squares, this can 
trigger a catastrophic failure of most of the network. Such extreme fragility is 
why the set of all frozen configurations has measure zero in infinite systems in 
spite of its non-zero entropy density. 

Before moving on to more complicated KA models, it is important to note 
that the whole proof of irreducibility is based on two ingredients: all config- 
urations that contain a special frame of holes belong to the same irreducible 
component; and the frame can be created by starting from a small nucleus of 
vacancies and expanding to larger sizes by satisfying at each step requirements 
that becomes less and less severe probabilistically. Therefore, the basic result 
for the simple square lattice model we have considered thus far is not strongly 
dependent on the exact set of dynamical constraints: in the next subsections 
we show how it can be adapted to all KA models on hypcrcubic lattices. 
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Figure 5: Construction of a large framed square from a partial set of framed 
subsquares: the two configurations can be connected by a sequence of moves 
similar to those of figure starting from the corners that connect two empty 
edges one can move lines of vacancies up, down, right and left to reach the final 
configuration. 



2.1.4 Exponential convergence of frameability of large systems 

We now show that the probability a configuration of a hypercube is frameable 
approaches one (at least) exponentially fast in its linear size L, i.e. 

HA, p {F)>l-Cexp(-^j (11) 

for L > E u (p). 

In addition to its direct relevance for the square lattice model, this result 
will be needed for analyzing other KA models. It is instructive to consider two 
ways of proving the exponential convergence. 

(i) The first is via a percolation-like argument, in the same spirit as that used 
in reference [42]. Consider initially a square lattice A of linear size L, 
divide it into a set of smaller squares of linear sizes I and focus on one 
of them, Ag. If the four neighboring squares of A; arc framed, then the 
31 x 31 square that includes all these small squares can also be framed (see 
figure |SJ); other combinations of the nine subsquares also suffice to make 
the larger square frameable. 

In this way one can grow large frameable squares. The crucial observation 
is that if the (possibly) unframeable subsquares do not span across the 
larger square, then the latter can definitely be framed. This condition is 
simply the lack of percolation of unframeable squares. Thus 

1 - HK h , p (F) < P perc {L/l^ Kup {T),2) (12) 

where P perc {l,r,2) the probability of conventional site percolation for a 
square lattice of linear size I, with occupation probability 1 — r. From 
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site-percolation estimates [42] 

P P erc(l,r,d) < f^l-r)^" 1 = (13) 
Therefore, if for the sub-squares 

/iA !>P (^) > i - Y e ( 14 ) 

equations 1)12)1 and 1)13)1 imply that for the large square 

/M.pCT 7 ) > l-Ce- L / 1 (15) 

with a positive constant C. Since 1)14)) holds for / = S — by the definition 
of the crossover length — the desired exponential bound, tf 1 1|) . is proven. 

(ii) The second method is via a renormalization group (RG) approach (see [60] 
for a similar argument in the context of bootstrap percolation). Consider 
a square lattice of linear size L and divide it into four squares of linear 
size L/2. If at least three of the four subsquares are frameable, then the 
large square is also frameable. This implies an iterative inequality: 

/M 21lP (^) > (pA L jT)f + 4 {^ L , P {F)f (1 - Ml A?)) ( 16 ) 

If 

1-»a l AF)< ^ (17) 

then on the next scale up 

VA 2L A^)> 1-10(1 -MA L ,p(^)) 2 (18) 

which implies that condition 1)17(1 holds also for 2L. We can thus iterate 
the inequality (fT£|) to obtain 

l*A*> L A?) > 1 - 1/10 (10(1 - /iA^p(^))) 2 " > 1 - ^e" 2 " (19) 

Since for L > S(p) condition (|17|l holds, equation 1)19(1 gives the desired 
bound (fT5|l for large square lattices. 4 

2.2 Irreducibility for d=3, s=2 

We now turn to the double- vacancy-assisted cubic lattice case, d = 3 s = 2, that 
was originally studied by Kob and Andersen, and prove that configurations are 
almost surely in the same large irreducible component in the infinite system 
limit for any density. This case can be analyzed by the general discussion of 
section l2~4l however the study of this more highly constrained three-dimensional 
case is useful as it introduces, in a readily visualized situation, the technique we 
use to extend ergodicity results from smaller to larger values of d and s. 

4 Note that we have only considered L = 2"S2,i(p) for integer n, but the results can be 
readily extended to all L. 
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Consider a cubic lattice A £ 1? with size L 3 . Define framed configurations 
as those with all the boundary sites empty and define frameable as those that, 
by an allowed sequence of elementary moves, can reach a framed configuration. 
Again, frameable configurations with the same number of particles belong to the 
same irreducible component. This can be checked by noticing that the bottom 
and top planes of vacancies can be raised and lowered performing the same 
"sandwich procedure" as for the d = 2, s = 1 case, but here using empty planes 
instead of empty rows. 

A lower bound on the probability that a configuration is frameable, ^A jP (^ r ), 
can be obtained by the following expansion argument. Consider a cube of linear 
size four, which has at its center an empty cube of linear size two. If adjacent to 
each face of the internal cube there is a four by four square that is frameable in 
the d = 2, s = 1, — single-vacancy-assisted planar — sense, the whole cube of 
size four is frameable in the desired three-dimensional double-vacancy-assisted 
sense. This procedure can be iterated to grow an L by L by L frameable con- 
figuration from an empty cube of size two, by requiring that in each subsequent 
shell all six of the square faces are frameable. 5 . 

Letting /iA,p(-? r0 ) be the probability of frameable configurations grown from 
a nucleus at the origin, a series of lower bounds for 0A,p (•?"") follows beginning 
from a completely empty X 3 nucleus for general X 

L/2-1 g 

MA,p(-n>MA, P (^°)>(i-p) x3 n (aLp^)) yx ( 2 °) 

l=X/2 

with A2L is a 2L by 2L square corresponding to one of the expansion faces. From 
the previous subsection we know that the frameable probability of a square, 
^''(f), is exponentially close to unity for L 3> £2.1(0). As in the case of 
bootstrap percolation [41] choosing X >> £2.1(0) in i|20[) makes the pre-factor 
dominate over the product and we conclude that 

for L S> £2,1(0); this yields 

6,2(0) ~ £2,1(0) (22) 

as the characteristic length scale above which the frameability conditions are 
met with high probability. 

By considering the O (L 3 /^ 3 ) possible initial nuclei in a large L 3 cube, we 
can show, by analogy with the two-dimensional case, that lim/^oo 0A,p(-7 r ) = 1 
thereby establishing an upper bound for the crossover length £3,2(0) above 
which the large irreducible component dominates: 

S^, 2 (0)^6,2(0)(l-0)- S ^ (p)3/3 . (23) 
Up to subdominant corrections, this means 

log£^ 2 (0)~[£^(0)] 3 (24) 

5 Actually we should require a little more than that all the square faces be frameable: one 
needs two frameable square of size L by L, two of size L by L + 2 and two of size L + 2 by 
L + 2. We skip these details whose effects are subdominant compared to those we compute. 
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2.3 Irreducibility for all d, s = 1 

We now extend the results derived in the previous two subsections. First we 
consider s = 1 on general hypcrcubic lattices. 

As above, define framed configurations to be those with all their one-dimensional 
hyperedges empty and define as frameable those reachable from framed config- 
urations. In the three-dimensional cubic case, for example, the relevant edges 
are the twelve edges of a cube. As in the cases already discussed, all framed 
configurations — and therefore also the corresponding frameable ones — belong 
to the same irreducible component. This can be seen, for example, in the cubic 
case, by considering permuting the occupation variables of two neighboring sites 
{i, j} which lie in a plane parallel to the "bottom" plane. By starting from the 
corners, the bottom frame can be raised till it frames the plane that contains 
the pair {i, j}. Then, after applying on this plane the same sandwich technique 
used for single-vacancy assisted motion on a square lattice, the desired permu- 
tation can be performed and the frame then reexpanded and returned to the 
bottom plane. 

The growing technique that enables the construction of larger frameable con- 
figurations starting from an empty nucleus, can be generalized via the following 
observation. Given a frameable hypercube of linear size L, to expand it to size 
L + 2, 2 d_1 vacancies are needed in the subsequent shell for each of the 2 d hy- 
perfaces. Therefore, the probability that a configuration is frameable, fJ^ p {^F), 
can be bounded from below by M/fp^ ); with 

00 / 2 d " 1 -l ]d _ 1 . \ 2 

where the product is over even I. In the limit p — > 1, L — > oo 

MS(^)~expf--^) ) (26) 

where 

c(d) = -2- V - 1) f log ( 1 - y ( 2 E 1 t] T l ) ) (27) 

As for the square case, there exists a length 

d, 1) « (- log(l - p)/(l - p)) 1/(d - 1} (28) 

such that ^^^(JF ) varies little with L for L ^> £. Therefore, by considering the 
O (L d /£ d ) possible distinct positions for the nucleus, an upper bound for the 
cross-over length follows 

ElM ~ exp ( C{d] ] ) (29) 

As explained at the end of section 12.11 requiring that the number of possi- 
ble positions for framing nuclei times the probability that the configuration is 
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frameable about a given nucleus is of order one, is not enough to conclude that 
HA.pi^ 7 ) = 1- Again, the conditions that a configuration is frameable about two 
different nuclei, are not independent events. However, provided the nuclei are 
sufficiently far apart, these events are almost independent since the frameability 
requirements become weaker and weaker for larger sizes. Indeed the requirement 
to have at least 2 d — 1 vacancies on each of the 2 d faces of a frameable hypcrcube 
only amounts to 2 d — 1 vacancies among L d_1 sites and thus, at fixed density, 
this becomes less and less restrictive at larger sizes. 

We outline an iterative argument by which the desired previous result can 
be proved (see [60] for a similar argument for bootstrap percolation). Divide the 
lattice into {Ljl2) d hypercubes of linear size I?,. At the center of each of these 
hypercubes consider a smaller hypercube Ai of size l\ = l%/2. The probability 
that the lattice is frameable can be bounded using the probability that the lattice 
can be made frameable starting from one of the (Ljl^f hypercubes of linear size 
l\. We thus have a bound ^a. p (JF) > P1P2P3, in terms of the probabilities of 
three events: Pi, the probability that at least one of the hypercubes of size h is 
frameable; P2, the probability that this frameable hypcrcube can be expanded 
until the size I2 (requiring 2 d ~ 1 vacancies for each shell from size 1% to size I2)', 
and P3, the probability that the frameable hypercube of size I2 can be expanded 
until the size L requiring2 d_1 vacancies on each shell from size I2 to L excluding 
the regions of space occupied by the other small hypercubes of size l\. The 
following results can be readily derived: 



Pi = ^A^)) 4 - !- ex P -MA l!P (^ 



L 



* = n 1 e (tk *(i-p) 4 

l=h \ i=0 

p 3 = f[(i- 2 j:\^r)p m)d ~ i -^-pY 

l=U \ 1=0 



2" 



(30) 




By choosing l\ m £(p, 1, d), these give 

Pi fn 1 — cxp ^— exp 
P 2 w 1 

P3 « 1 (31) 

so that /LiA jP (J r ) is of order unity as long as exp(— c(d)/(l — p) 1 /^ -1 ^) (L/l 2 ) d ^> 
1, i.e. for L > 3^(p) with ^lip) S iven b y Eo.^25j>. 



2.4 Irreducibility for all d with s < d 

We now consider KA models with general parameter < s < d on a hypercubic 
lattices of dimension d. These are the non-trivial cases. [As mentioned earlier, 
or s > d, any completely occupied d-dimensional hypercubes are frozen so that 



21 



these KA models on hypercubic lattices are non-ergodic for any density p > 0. 
At the opposite extreme, s = is the normal lattice gas which is trivially 
ergodic] 

The s = 1 case was considered in the previous subsection. We now show how 
the general case can be analyzed, and irreducibility proved, via a generalization 
of the iterative procedure used in section |2~2*1 to extend the d = 2, s = 1 results 
to d = 3, s = 2. 

Define framed configurations of (i-dimcnsional hypcrcubes as those with all 
hyperedges of dimension s empty. [More formally the configurations with no 
particles on the hyperplancs: 

{x El =1,L}X {x E2 =1,L} X ... {x Ed _ s = l,L} (32) 

where Ei, . . . Ed- S for all sets of d— s of the indices 1, . . . , d.\ Inside a frameable 
hypercube, any move of a particle to a neighboring empty site can be achieved 
by the generalization of the sandwich technique discussed previously. 

Consider a frameable hypercube of size /. If adjacent to each hyperface there 
is a d — 1-dimensional hypercube that is frameable in the d — 1, s — 1 sense, then 
the enclosing d-dimensional hypercube of size I + 2 is frameable. Therefore, to 
bound from below the frameable probability, fJ,A,p{^ d ' s ^), we need to estimate 
the probability that a hypercube of linear size X centered at the origin is empty 
and that expanding from it any subsequent shell is frameable in the d— 1, s — 1 
sense. This yields 

L 

MA*,p(-?fd, a )) > 0--p) X II M/wO^d-M-i)) (33) 

l=X + l 

where A^ _1 is a d— 1-dimensional lattice of linear size £. The argument given is 
section l2\l.4l for d = 2, s = 1 can be extended to the get that the convergence of 
the frameable probability to one is exponential in the linear size of the system 
and determined by s _ 1 (p) 6 

MA,p(^ (( i-M-i)) > 1 - Ccxp [L/SH-l-i] fOT L » SJ-Ld-itp)- (34) 

Hence, by choosing X S> d _ 1 (p), the product on i in is close to unity. 
From generalization of the arguments given for the s = 1 square lattice case, 
a non-zero probability of frameability around a given nucleus that converges 
rapidly to its asymptotic limit for large sizes implies a probability close to unity 
for frameability of sufficiently large hypercubes around some nucleus. More 
precisely that PK.p{^F(d,s)) approaches unity in the thermodynamic limit and is 
close to one for L>SJ s (p). Since the dominant small factor in (|33fl is the large 
power of 1 — p, and to make frameability likely one needs a system large enough 

6 Note that the argument for exponential convergence requires the a priori knowledge of 
convergence to unity for the probability of the frameble component. This has been proven 
for any d, s = 1 in previous section. Therefore 1341 holds for any d, s = 2 and the argument 
below leads to fiac.p = 1 for these cases. Therefore, one can apply again arguments in section 
12. 1.41 and obtain exponential convergence for all the cases s = 2, i.e. 1341 for s = 3. And so 
on iteratively. 
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to contain at least one frameable hypercube of size greater than s l , an 
upper bound on the crossover length is 

[(l-p)- 3 "^"- 1 '- 1 )*]* . (35) 

From the earlier results for s = 1 in general d > 2, wc know l E^ 1 (p): Eq. (|29|l . 
Therefore, by induction, equations (|35fl yield upper bounds for the crossover 
length for general s, d with 1 < s < d — 1. The density dependence of this 
upper bound, S", is hence 

^(p,s,d) = e X p°s {1 °^jl_ g) (36) 

where exp os denotes the exponential function iterated s times. Note that the 
factors of ln(l — p) that appear in 1351 can be taken into account, as far as the 
upper bound, by a slight modification of C(d, s). 

2.5 Ergodicity for all d with s < d 

Thus far, we have only considered the irreducibility of configurations in infinite 
systems. In this section we prove that the almost-sure existence of an irreducible 
component in the thermodynamic limit implies ergodicity (i.e., the third step of 
the procedure outlined in section l2~T|) . We follow the same strategy as in section 
4 of [43]. Since the proof is the same for any of the models, we temporarily drop 
the indices d and s. It is first necessary to introduce some mathematical notation 
and give a rigorous definition of ergodicity for infinite systems. 

Define Q\ = {0, 1}' A ' as the configuration space of the hypercubic lattice 
A G 7L d and configurations rj as the elements of fi, namely sets of occupation 
numbers rj x G {0, 1} for x G A. The dynamics of KA models is a continuous 
time Markov process with generator C acting on local functions / : Ha — > K as 

£f(v)= E Vx(l-Vy)c xy (v)[f{V xy )-m} (37) 

{x,y}cA 
\m,y\=l 

with the sum running over nearest neighbor pairs, x, y. The quantity in square 
brackets, 

V xy f = f{rf v ) - f(v) , (38) 

acts to change the occupation of sites x and y so that the occupation number of 
the zth site in configuration rf v is related to the occupations in the configuration 
X] by 

{r\ y if z = x 
if z = y (39) 
r] z if z ^ x, y ; 

The jump rates c x . y (r]) encode the constraints imposed by the KA rules, namely 
, . _ J 1 if nj (17) < m and n^(r?) < m ( 

Cxy[V) — \ otherwise l 4U i 
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where 



n 



(41) 



zSA, z=iy 
|o?— ar|=l 



i.e. n^.(rj) is the number of occupied neighbors of x, excluding y. 

The dynamics preserves the number of particles, [i.e. the subspaces with 
fixed number, TV, of particles £Ia,n := {rj G Qa Y^xeA^ = ^} are invariant 
and the uniform measure va,N on fixed N subspaces is invariant]. However, 
due to the vanishing of certain rates for m < 2d — 1 (s > 1), the configuration 
space, r^A.TV, of finite lattices breaks into disconnected components: the Markov 
chain is hence reducible so that the process is not ergodic on these subspaces. 
On the other hand, for the infinite lattice A = 1 d . the process satisfies de- 
tailed balance with respect to the trivial Bernoulli product measure, /j, p , at any 
density p. Moreover, from the previous section, we know that there exists an 
irreducible component — the set of frameable configurations — that covers the 
configuration space. However this is not enough to conclude that the system is 
ergodic. Indeed, while for finite state systems irreducibility of the Markov chain 
implies that there exists a unique invariant measure and the system is ergodic 
with respect to it, this is not a priori true in infinite systems. For example, Ising 
models below their critical temperature are irreducible but not ergodic in zero 
field because of the coexistence of the two phases each of which lasts for infinite 
time. To establish crgodicity for infinite systems, one needs to prove [44] that 
the long time limit of all correlation functions — weighted averages over the 
probability distribution at time t, Vt = Vt ~ exp Ct — approaches those of the 
Bernoulli product measure for almost all initial conditions, rj : more precisely 
that 



where / G L2([i p ) if / dfi p (r])f 2 (r]) < oo and n P (f) is the equilibrium average of 
function /, namely 



and VtfiTj 1 ) is the correlation function f(rj) averaged over all possible histories 
up to time t starting with initial configuration r] 1 . 

By the spectral theorem, convergence to the equilibrium measure occurs if 
and only if zero is a simple eigenvalue of the generator of the dynamics, £, i.e. if 
the only functions, /o, in L2 (/•*/>) f°r which Cfo(rj) = are constant on almost all 
configurations, i.e. on all except possibly a set of measure zero [44]. Note that if 
the system can be trapped for an infinite time in some region of the configuration 
space, this would not be true: the characteristic function of this region would 
be a non-constant eigenvector of C with zero eigenvalue. This corresponds to 
the natural idea that ergodicity breaking is related to the existence of regions 
of the configuration space that are effectively disconnected. In the Ising model 
case, the function f(rj) that is equal to unity on configuration 77 if a majority of 
the sites have up spins, and equal to zero otherwise, is such an example. 

The strategy of the proof of crgodicity of KA models on infinite hypercubic 
lattices Z d with parameter s such that s < d is the following. (1) First we prove 





(42) 
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that, if fo is eigenvector of C with zero eigenvalue, fo{v) = fo(v x ' v ) almost surely 
with respect to \i p for any pair of sites {x,y}, i.e. ^ P ([f$ v — fo} 2 ) = where 
fo V {r]) = fo(v x ' y )- This part of the proof uses as a key ingredient the existence 
of an irreducible component which has unit probability in the thermodynamic 
limit. (2) Then we use De Finetti's theorem [45,46] on the decomposition of 
symmetric measures on product measures to conclude that any such /o is in fact 
constant. 

Let us prove (1). By enumerating framcablc configurations, t] , rj 2 , ... rj n , . . . , 
we can rewrite the event that a configuration is frameable as T := U n >i 7 ?" arL d 
define the function Ijr to be one if rj 6 T and zero otherwise. The result 
fJ- P (J~) = 1 established in previous sections implies 

Miv^/] 2 ) = M[v^/] 2 M < E M[v*„/] 2 V) v/ e L 2 ( Mp ) (43) 

n=l 

where the exchange operator V xy was defined in |JSSJ. From the properties of 
frameable configurations, for any configuration rj G T one can always find a path 
r?, rj 1 , • • • , ?7 with ry M = rf v that connects ?7 to rf y through allowed moves, 
i.e. with rf +1 = rf'' zw for some pair of neighboring sites {z,w} and c zw (rf) = 1. 
By telescoping sums and the Cauchy-Schwartz inequality 7 , each term in the 
sum (|43[) can now be bounded from above by a sum of elementary exchanges 
allowed by KA rules, i.e. terms jJ. p (c zw [\7 zw f] 2 ). Note that we can introduce 
the jump rate in the upper bound only since, as written above, we can choose 
elementary moves with rate equal to one (never to zero). This is thanks to the 
fact that the frameable component has unit probability, thus we could restrict 
the mean in (|43|l as a mean over the frameable component, where we know 
that a path of allowed moves exists for sure. The hypothesis of sentence (1), 
namely that there is some fo that is an eigenvector of C with zero eigenvalue, 
implies that fi p (foCfo) is zero. Moreover, the symmetry c xy {rf) = c yx {rj) — the 
microscopic reversibility — implies 

H P {foCfo) = -\ J2 VpMl)lV xy fo} 2 ) ■ (44) 

Since the right hand side of l|43|) can be rewritten as a sum of such terms, the 
fact that ji p (foCfo) = implies jJ, p ([V x y fo] 2 ) = and completes the proof of 
(I)- 

We now turn to the proof of (2). Consider the measure Ji*° defined as 
fjf°(r]) = /q (ri)(J. p (r])/ ji p (fQ ). Since ji p is symmetric under exchanges of particles 
and vacancies, i.e. fJ, p (r]) = ji p (rj x - v ), and fo{"n xy ) = foiv) on almost all the 
configuration space, fi^° is also symmetric, i.e. ^"{rj) = jl^°(r] x ' v ) for any r] 
and any pair {x,y}. Therefore, by the De Finetti theorem [45,46], it can be 

7 Note that for any frameable configuration the existence of a path that joins rj and rj x ' v is 
guaranteed. However, its length can be arbitrarily large on some rare configurations. There- 
fore, strictly speaking, one cannot directly apply Cauchy-Schwartz inequality which would 
pull out a non-bounded factor proportional to the length of the maximal path. However, by 
dividing framcablc configuration on sets characterized by the length i of the minimal path (to 
perform such move), it is possible to bound each term in the r.h.s. of 1431 as an infinite sum 
(on t) of terms of the kind i multiplied by terms pi p (c zw [S/ zw f] 2 ). This, together with the 
observation (see below) that each of this terms must be equal to zero, is sufficient to conclude 
that the l.h.s. in 1431 is equal to zero. 
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decomposed on product measures corresponding to different densities. However 
measures corresponding to different densities are concentrated on different sets 
of configurations and, hence, all the coefficients of the decomposition have to 
be zero except the one corresponding to p, p therefore /o (^V/vC/o) = 1 which 
implies that fo (tj) is constant 8 and equal to the coefficient of the decomposition 
corresponding to p, p 9 . 

3 Quantitative results 

In previous sections we have proven that for any parameters d and s, with s < d, 
crgodicity holds in the thermodynamic limit. Moreover, we have determined an 
upper bound H^ s (p) for the density dependent cross over length, Sd iS (p), that 
separates two different regimes on finite volume hypercubes A £ Z d of linear 
size L: if L — > oo and p — ► 1 with L > 5 Sj£ ;(/9), than p\^ p (M) — > 1; if L — ► oo 
and p — ► I with i > S S ( j(/o), than /za, p (.M) — > 0, where .M is the maximal 
irreducible component (i.e. set of configurations on A connectable by possible 
path which has larger probability w.r.t. Bernoulli measure at density p, pa.p)- 
In other words, for L <C S Sj( j(p) the configuration space is dominated by a single 
crgodic component, while for L <C S Sj( i(p) different components contribute. 

In the next subsection we find a lower bound S l d s (p) , which has the same 
form of density dependence as the upper bound. In the subsections 13.213.31 
we show that for specific cases, the exact asymptotic behavior of the crossover 
length can be obtained — including the numerical coefficients. 

3.1 Bootstrap percolation and lower bounds on crossover 
length 

In this subsection we explain how one can obtain a lower bound on the crossover 
length, S, thanks to a relation among KA models and the well studied problem 
of bootstrap percolation. 

Consider a random configuration at density p on a d-dimcnsional hypcrcubic 
lattice A. One can isolate frozen subsets (possibly empty) of such a configuration 
by the following deterministic procedure. First remove all the particles that have 
no more than m neighbors, then iterate this procedure until no more particles 
can be removed. Indeed, as already noticed by Kob and Andersen for the d = 3 
s = 2 case, all the particles (if any) that remain at the end of such procedure 
are frozen: namely they can never move under KA dynamics starting from 
the original configuration 10 . This procedure is simply related to conventional 
bootstrap percolation ( [47], [48], [49]) which corresponds to starting from a 
configuration at density p, itcratively adding particles in empty sites that have 

8 Since this procedure can be repeated for any linear combination of fo(v) with the constant 
function, we obtain that the square of any linear combination of foiv) with the constant 
function is a constant which implies that fo(r]) is indeed itself constant. 

9 Note that for Ising models at low temperature in zero field we cannot invoke such a unique 
decomposition (the hypothesis for De Finetti's theorem do not hold since the equilibrium 
measure is not simply a product) and therefore ^ p {[}q ,v — fo] 2 ) does not imply that /o is 
constant. 

10 Note that the converse is not guaranteed: a priori there could exist sets of particles that 
are thrown away during the deterministic procedure, but they are forever blocked under KA 
dynamics. 
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fewer than m neighbors, and considering any clusters of empty sites that may 
remain at the end of this process. In other words, by exchanging particles and 
vacancies, and hence p with 1 — p, one recovers the usual bootstrap procedure. 
Let p B be the probability that, starting from a random configuration on an 
hypercube A 6 Z d of linear size L, a cluster of particles remains at the end 
of the modified bootstrap procedure we are interested in. The rigorous results 
in [42] for usual bootstrap percolation yield the existence of a crossover length 
E B (p,d,s) = exp os (K(d,s)/(l — p) d ~ s ) such that if one sends simultaneously 
p — * 1 and L — ► oo with L growing faster (slower) than 5 s , p B p goes to one 
(zero). Furthermore, the convergence to unity in the regime for L > S B (p, d, s) 
is exponentially fast. 

To our knowledge only upper and lower bounds on the constant K(d, s) are 
available in the general case, while the exact value has been recently determined 
[50] for the case corresponding to our d = 2, s = 1. For the non-trivial models, 
s < d, if clusters of particles remain, they must be system-spanning clusters (i.e. 
clusters that connect two boundaries of the system): in infinite systems, such 
frozen clusters must be infinite. Since on any finite lattice two configurations 
with different frozen sets belong to different irreducible components, the likely 
presence of such system-spanning clusters of frozen particles implies immediately 
that when p — » 1 and L — ► oo with L increasing slower than S B (p, d, s), the 
probability of the maximal irreducible component for KA model (p\ tP (M) is 
strictly smaller than one. This corresponds to the following lower bound for the 
KA crossover length: 

S d , s (p)>S^(p) = S B (p,d, S ) . (45) 

Note that the form of the density dependence of this lower bound (|45J) is the 
same as that of the upper bound, ()36|l . 

As mentioned above, for the case d = 2 s = 1 recent results have been found 
in [50] with a sharp value of the constant: K(2, 1) = 7r 2 /18. On the other hand, 
from equation flJJb we find (7(2,1) = - J* dy ^s{i-v+y^v) _ 4 4; thcrc f orc 
K(2, 1) < C(2, 1) and Sj; i(p) > Z l 2 \{p)- In the next subsection we show how, 
from a different framing technique, it is possible to determine a stronger upper 
bound for S2,i(p) which turns out to be equal to the lower bound from bootstrap 
percolation. 

3.2 Optimal framing and exact crossover length for d = 

2, s = 1 

We now show how exact results for the crossover length of the single-vacancy 
assisted (s = 1) KA model on a square lattice can be obtained by an optimal 
framing construction. 

Define a W x H rectangle (with W + H even) to be optimally framed if it 
has \{W + H) + 1 vacancies arranged with one in a corner and the others on 
alternate sites of the two perpendicular sides that intersect at that corner, plus 
any number of additional vacancies, see figure Optimally frameable configu- 
rations are those that can be reached from an optimally framed configuration 
with allowed moves. 

One can show that all optimally framed configurations of a rectangle with 
the same particle number belong to the same irreducible component. Hence all 
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Figure 6: A 7 by 7 optimally framed configuration 
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Figure 7: Basic moves by which the three vacancy corner can move along the 
alternating-vacancy row of an optimally frameable rectangle. 



the optimally frameable configurations belong to the same irreducible compo- 
nent. This follows by showing that within an optimally frameable rectangle any 
nearest neighbor pair of particles can be exchanged and all the other particles 
returned to their starting positions. The sequence of moves to perform a generic 
such exchange can be constructed by considering the basic moves in figure 
which allow one to move the nucleus of three vacancies along the row containing 
alternating vacancies, and those in figure |H1 which enable one to lower and raise 
the row of alternating vacancies through the lattice. 

It is straightforward to check that with fewer than 1 + ^(W + H) vacancies, 
large scale rearrangements of the particles in a rectangle are not possible: this 
suggests that this framing is, at least in some senses, optimal. 

Following the same strategy as for the simple framing used previously, we 
now show how larger optimally frameable rectangles can be constructed. Con- 
sider an optimally framed W x H rectangular region embedded in a larger 
system. The following statements can be checked by direct inspection: if there 
is a vacancy in a line next-nearest neighbor to one of the rectangular edges 
parallel to either of the directions x or y, the rectangle can be expanded to a 
W x (H + 2) or (W + 2) x H, respectively, optimally framed rectangle; if there is 
a vacancy in the line segment next to one of the rectangle's edges it is likewise 
expandable into a W x (H + 1) or (W + 1) x H framed rectangle; and if there 
is a vacancy next nearest neighbor along a diagonal from one of its corners, the 
rectangle is expandable to a (W + 1) x [H + 1) optimally framed rectangle. 

Starting from a nucleus of three vacancies in a two-by-two square, the de- 
scribed expansion procedure can be iterated to grow larger frameable rectangles 
as long as the needed vacancies are present at each step. We must now estimate 
the probability that all the needed vacancies to construct an infinite frameable 
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Figure 8: Basic moves by which the row of alternating vacancies in an optimally 
frameable rectangle can be lowered. 
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region centered on a chosen site are present. Define Q e k as the probability that, 
given a k — 2 x I optimally frameable rectangle, it can be expanded to the right 
to an optimally frameable k x I rectangle which now includes the (k — l)th and 
fcth columns. From the above observations, we obtain the following recursion 
relations 



lk+2 



-vi\ 



f k e-^(l e 



(46) 



where we have approximated, in the high density limit of interest, p e ~ e e ^ p ) 
(and used the notation v = 1 — p). Defining the ratios of successive Qs by 



R 



i _ 



Qk+l 



(47) 



we obtain 



Rl +1 = (l-e- vi ) + 



-v£\ 



R 



(48) 



Since the least unlikely way for the rectangle to expand is roughly isotropically 
— as can be shown — , we expect k sa I at large scales . Moreover, for k and 
I large the ratios of probabilities for successive expansions along the same axis 



will vary slowly, therefore a natural - 

n k+l ~ 

By substituting (gSJ in we find 



- and checkable — ansatz is 
R[ « R{i) 



R{1) 



where 



£ 



-vl 



(49) 



(50) 



(51) 



As the rectangle has to be expanded to infinity in all four directions concurrently, 
this yields an estimate for the probability /^^(O) that a specific nucleus of 
three vacancies can be expanded to an infinite optimally frameable rectangle 



exp 



(52) 



which, by replacing the sum over £ by an integral and changing variables, yields 



2c n 



AWp(O) ~ ex P~Y3^ 



with the constant given by 



d£ 



in 



IS 



0.55 



(53) 



(54) 



The knowledge of Poo. p{0) gives, as in previous sections, an upper bound for 
the crossover length 



-2,1 (p) < exp 



(1-P) 



(55) 
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Figure 9: Rectangular unit cell for the basic planar structure for the d = 3, 
s = 2 case 

Since the constant Coo Ij54(l has exactly the same value as the bootstrap constant 
K(l, 2) in [50], the upper bound for 3 obtained by the above framing procedure 
coincides with the lower bound provided by bootstrap results, thus giving the 
exact asymptotic high density behavior of the crossover length, up to subdomi- 
nant pre- factors. This proves that the optimal framing procedure captures the 
dominant mechanism that restores ergodicity in large systems at high densities; 
this justifies referring to it as optimal framing. 

3.3 Optimal framing for d=3 s=2 

We have shown that the framing procedure devised in section \2. II for the two- 
dimensional case is not optimal: indeed it does not capture the dominant er- 
godicity restoring mechanisms and therefore gives a non optimal bound for the 
crossover length, 2. (see also section QJ. However, we were able to construct a 
different framing procedure which does capture the dominant mechanism. This 
was done through the construction of a lower-dimensional basic structure (lines 
of alternated particles and vacancies) in which each particle is blocked but only 
barely so: an additional lowcr-dimcnsional vacancy cluster (an extra vacancy in 
the square lattice case) can be moved throughout the whole structure. 

We now construct a framing that we conjecture is optimal for the d = 3 
s = 2 case. Consider a cube of linear size L which is completely filled with the 
exception of a plane L x L on which the configuration is given by the repetition 
of the unit cell of eight sites drawn in figure|5| It is possible to check that all the 
particles on such a plane can never move. However, if one adds two vacancies 
in place of e.g. the first particles on the two rows of one of the cells, then all the 
particles on the plane can be moved. Consider a cube which has three orthogonal 
faces with the above defined basic planar structure. We conjecture that with an 
extra set of O(L) vacancies associated with the linking edges and 0(1) with the 
corner, the planar structures can be moved all across the cubic lattice and allow 
any exchange inside, i.e. they play the same role as the completely empty faces 
of frameable cubes defined in !2.2l If this is correct, by requiring 9/8L 2 + 0(L) 
vacancies one can construct a configuration for a cube of linear size L in which 
all exchanges are allowed. We conjecture this to be the minimal requirement, 
i.e. the above defined framing procedure to be the optimal one. 

Note that the calculation of the probability for such a framing procedure 
to be expandable should yield an upper bound on the crossover length which 
is better than, although of the same form as, the one we obtained in (|23Jl . 
since fewer vacancies are required in successive layers. However, a rigorous 
lower bound — and hence the putative exact value of C(3, 2) if the framing 
described here is indeed optimal — is not known; the exact value of K(3, 2) for 
the corresponding bootstrap percolation has not been determined. 11 

11 Note that knowledge of C(3, 2) from framing arguments would give a new upper bound on 
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4 Physical picture of the cooperative KA dy- 
namics at high density 

In this section, for completeness, we summarize the physical picture that we 
derived for the high density dynamics of the KA. This has already been discussed 
in [37] and the derivation of some results on the dynamics will be presented in 
a forthcoming publication [51]. We focus for simplicity on the s = 1 d = 2 case 
leaving discussion of the natural generalization to any of the s < d KA models 
on hypercubic lattices to the end. 

As already noted, in the interesting cases there are no finite clusters of 
vacancies that can freely move in an otherwise completely filled lattice. However, 
in a frameable region of size I — £ ~ ln(l/(l — p))/(l — p), we know that: (1) all 
particle configurations inside the region can be reached from each other via an 
allowed (although long) path of elementary moves; (2) a frameable region can 
move one step in any direction if it can expand one step in that direction (it 
expands and then leaves a vacancy behind on the opposite side when it moves). 
(3) a frameable region of size I = £ ~ ln(l/(l — p))/(l — p) is likely to find at 
least one vacancy in each of I consecutive line segments along each of the lattice 
directions. Thus frameable regions of size £ can expand and, hence move, £ 
steps, a distance similar to their diameter. For this reason, we refer to them as 
mobile cores [37], although it may be a bit misleading, they have been called 
defects in other contexts. 

The motion of frameable cores is, on long scales, like simple diffusion in a 
random medium. The size of mobile cores, £ = £ ~ ln(l/(l — p))/(l — p), is such 
that the random environment seen by the cores is the minimal one in which 
diffusion is possible. For smaller values of £ the probability to find the neces- 
sary vacancies in the neighborhood of a core is too low to guarantee that it can 
diffuse. In contrast larger frameable regions with £ > £ could diffuse, but these 
would contribute less to the long-time dynamics because they are both rarer 
and move more slowly than cores of size £. [Indeed, larger frameable regions 
will tend to decompose into one (or more) minimal mobile cores, and leftover 
vacancies that cannot move without assistance of a mobile core.] 
The diffusion coefficient of a tagged particle, D$, one of the key physical quan- 
tities analyzed in simulations and experiments, will be given by the spatial 
density of the cores divided by the typical time scale for their motion, assum- 
ing approximatively independent core diffusion. The density of cores of size 
£ ~ ln(l/(l — p))/(l — p), is just the frameable probability around a given ori- 
gin, /.tA, p (^ 7 ) ~ The form of the asymptotic behavior of S in the high 
density limit was found for any s < d, and, for the d = 2 s = 1 case, the coeffi- 
cient that enters In S was calculated exactly see section [3] 

The typical timescale on which cores move has been discussed in [37]. Since 
all frameable configurations are equiprobable, this is proportional to exp(AS'), 
the ratio of the number of accessible (frameable) configurations of the core to 
the number in the most severe bottleneck in the configuration space of the 
core [37,51]. The worst case scenario for a minimally frameable region of size 
£ would occur if the bottleneck corresponded to a single configuration, (the 

the value of K(3, 2) for the Bootstrap problem and perhaps a hint of its exact value (recall that 
K(2, 1) coincides with the constant coming from the correspondent KA problem calculated 
with the optimal framing). 
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framed one) leading to AS = Stotal — hi£! since the total number of frameable 
configuration is of order £1. However, as explained in [37] (and to be detailed 
in [51]) the bottleneck is not nearly as tight: cores do not have to go through 
framed configurations in order to rearrange or to move. Numerically we have 
found that the time scale increases more slowly than exponentially in £ and 
its growth is compatible with an analytic argument for the entropy bottleneck, 
which yields at least a lower bound. This result is obtained by noting that, in 
order to fully equilibrate, an optimally frameable square of size I has to pass 
through configurations with the nucleus (the three vacancy element, see Fig. El 
and section I3.2J1 in a comer of the square. Using a transfer matrix technique 
we obtain for I x £ optimally frameable squares an entropy difference between 
configurations with the nucleus in the center (which are typical) and those with 
the nucleus in the corner asymptotically equal to AS w Tv* + a\n£ + C with 
T = 2\/ 6 + \/22— 2\/3 = 3.075 and a computable in principle. [51] This leads, 
via the observation that the worst blockages are typically when I ~ 1/(1 — p) 
to a typical timcscale for movement of optimally frameable regions 

lnrc^l/v/T - ~p (56) 

for the s = 1 square lattice model. We conclude that the dominant contribution 
to self diffusion of a tagged particle, D$, arises from the low density of mobile 
cores rather than the long time needed for them to move of the order of their 
diameter. 

The arguments given above can be generalized to any KA models on hypercu- 
bic lattices with s < d. The conclusion is that, to leading order in (1 — p)~( d ~ s \ 
the s— iterated logarithm of the self diffusion coefficient scales as the s— iterated 
logarithm of p^ p {T). Correspondingly, we expect the typical relaxation times 
t to scale with one over the density of cores, i.e. as the inverse of D$ 12 . In 
particular, thanks to the exact result in section l3~2l for the d = 2 s = 1 model 
we obtain 

lim(l -p)]nD s = tt 2 /9 (57) 

a prediction which we have successfully verified by numerical simulations [37]. 
On the other hand, for the case d = 3, s = 2, the dependence on density for 
different relaxation times r have been considered in [14, 24] and our predic- 
tion have been successfully checked, namely data are well fitted by the scaling 
t oc expexpC/(l — p) with C constant. 

Note however that there may be very substantial corrections to the leading 
asymptotic behaviour at obtainable densities. Indeed, it is known that there are 
large corrections to the asymptotic law derived for bootstrap percolation [47]. 
Therefore, it would be useful to compute sub-leading corrections to all the 
analytical expressions derived here. First, one should take into account the 
timescale for movement of a core, discussed above. For example in the d = 2 
s = 1 case corrections to (1 — p) log Ds from this will be large — - probably 
of order of y/1 — p. Furthermore, one should take into account corrections to 
the asymptotic form of the optimal framing probability: we expect these to 

12 Strictly speaking our results applies directly to Dg. The structural relaxation time could, 
and very probably will, affected by the slowest regions in the systems which instead are not 
relevant for the self-diffusion. Hence, we expect a structural relaxation time larger than l/Dg, 
i.e. a decoupling between structural relaxation time and self-diffusion as indeed the numerical 
results of Kob and Andersen already suggested originally [11]. 
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be negative as occurs for bootstrap percolation where the numerical results for 
Coo [47] correspond to roughly one half of the exact result [50]. [Note that re- 
cent progress on the calculation of such corrections for bootstrap percolation 
has been made.] [52]). A rough cancellation between such a negative correction 
to the density of the mobile cores, and the positive one arising from the typical 
timescale, might explain why, in the measured range of densities, In D$ seems 
to be closer to its asymptotic value than is the crossover scale S. 

The mechanism that we have identified for the long-time scale dynamics, 
certainly provides a lower bound on D$ for p — > 1 ; this could be probably be 
proved along similar lines [51]. But how do we know that there is not a faster 
mechanism which dominates? There are two reasons to strongly believe that 
there are no faster mechanisms: First, we have proven that the finite size cross- 
over between irreducible and non irreducible dynamics takes place exactly when 
the probability of having at least one core in the lattice crosses from small to 
large. This strongly suggests that a single core is necessary, as well as sufficient, 
to relax the system. This is exactly analogous to the behavior for a standard 
lattice gas in which D$ oc (1 — p) at high density and the change from frozen to 
unfrozen takes place when there is at least one vacancy in the region. Secondly, 
rigorous analysis of non-cooperative KA models [38] for which finite clusters of 
q vacancies can move freely, (such as s = 1 on a triangular lattice, for which 
q = 2) show that the self-diffusion coefficient scales as D$ oc (1 — p) q in the high 
density limit. The cooperative cases can be understood roughly in terms of a 
density dependent q: at any fixed density there are no free moving clusters but 
there is a characteristic sized group of vacancies — a mobile core — that can 
indeed diffuse since it typically finds enough outer vacancies in each direction 
(Of course the environment is random and there are some places in which a 
core cannot enter without changing its size. However these regions would not 
be important for the self-diffusion coefficient even if they were forever frozen, 
as happens for diffusion in a disconnected random environment when there is a 
percolating cluster.). At each density, one needs to have q oc £(p), in order to 
be sure that the diffusion of cores is possible: this yields a mobile core density 
of the form we have found. 

Up to now we have focused on the behavior at high densities, in particular 
the asymptotic value of relevant quantities in the limit p — > 1. However, as 
discussed in more detail in [37,51], for some lattices, especially with s > 1, 
the asymptotic behaviour will not be relevant except on extremely — perhaps 
inaccessibly — long time scales. In practice, much of the slowing down as 
the density increases could occur before this asymptotic regime. In particular, 
rapid crossover could occur near an almost transition — sometimes called an 
avoided transition [53] — near an almost-critical density ( which would depend 
on the lattice and the value of s). Such behavior could account, for example, 
for the apparent power law "vanishing" of the self diffusion coefficient observed 
as p — > p ~ 0.881 in simulations of the three-dimensional case with s = 2. 

A simple argument suggests that such a cross-over might exist. 

We again first focus on the two-dimensional case with s = 1. As already 
mentioned in [37] , in this case a small (fixed-size) cluster of vacancies can freely 
diffuse on a second-neighbor-percolating clusters of other vacancies. Therefore, 
when the vacancy density is above the correspondent percolating threshold, 
diffusion can occur without substantial coopcrativity. But the contribution to 
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the self diffusion coefficient from this mechanism shrinks to zero as the vacancy 
percolation threshold is approached from above (from below, in terms of particle 
density). At higher densities, the mechanism behind diffusion crossesover to the 
cooperative one discussed in this paper. 

For the three-dimensional case with s = 2, from the construction in section 13.31 
we expect that the approximate transition seen in numerics may well be related 
to percolation of surfaces on which the vacancies are nearly connected like those 
on the surfaces of optimally framed cubes. 

The above percolation arguments are only qualitative and we do not expect 
that the density at which the percolation transition takes place will correspond 
very well with the density at which the apparent transition takes place. In next 
section we present a quantitative analysis of an actual transition in KA models. 
We show that in finite dimensions, apparent transitions with rapid crossovers 
could arise as "ghosts" of actual dynamical transitions that occur on Bcthe 
lattices — or, cquivalently, within the Bethe approximation for real lattices. 
This suggests a way to study "avoided transitions" : by using the Bcthc lattice 
results and perturbing around this "mean-field" limit. 

5 KA models on Bethe lattices 

Bethe lattices are defined as a infinite loop-free graphs with fixed connectivity 
z = k + 1. The main motivation for studying statistical mechanical models 
on Bethe lattices is that such "lattices" are often considered to be a good ap- 
proximation of more realistic lattices, e.g. hypercubic, in the limit of high 
dimensionality: the Bethe lattice models are then considered to be a type of 
"mean field" — or more accurately "uncorrelatcd field" — approximation. But 
their tree-like structure enables analytic calculations via recursive procedures, 
As we will explain in detail, the behavior of KA models on Bcthe lattices 
is qualitatively different from that on hypercubic lattices obtained in previous 
sections: on Bethe lattices an ergodic/non-ergodic transition takes place at a 
non-trivial critical density. Nevertheless, exact results on these lattices, when 
combined with those for hypercubic lattices, should be useful both for better 
understanding of the latter, and as background for considering how mean-field- 
like scenarios for glass transitions might — or might not — be extended to real 
glasses, or at least to some more realistic models. 

5.1 Bethe lattices 

There are subtleties in the definition of Bcthe lattices that are important to 
note. Bethe lattices with connectivity k + 1 arc locally identical to the infinite 
size limit of random "c-regular" graphs with c = k + 1 [54] which arc uniformly 
drawn from the set of graphs with connectivity c for each site with neither 
multiple edges (no two edges joining the same pair of sites) nor loop-edges 
that join a site to itself. Since in the limit of large number of sites, with high 
probability there are no finite loops that go through a given set of vertices, a 
typical such random graph looks locally (i.e.: on any finite length scale) like 
a Cayley tree with a fixed branching ratio. [Recall that a Cayley tree with 
connectivity k + 1 is constructed by taking k + 1 rooted trees with branching 
ratio k and connecting all the k + 1 roots to a new site, which can be considered 
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as the origin.] Nevertheless, for macroscopic properties the presence of loops 
that exist in even very large random graphs is crucial since it induces geometric 
frustration and ensures a statistically homogeneous structure which prevents 
the pathologies that arise for Cayley trees for which a positive fraction of sites 
are on the boundary which causes extreme sensitivity to boundary conditions. 
A Bethe lattice corresponds to considering a Cayley tree of size L, focusing on 
a core of size I around the origin and taking the limit L — > oo before the limit 
I — > oo, or with the simultaneous limit with both lengths becoming infinite 
but -J ~~ * 0. Note that this is equivalent for models with local dynamics to 
considering first the limit L to infinity with random initial conditions, and then 
considering the long time limit: as such, it is quite physical. 
The following analysis holds both for random c-regular graphs with c = k + 1 
and for Bethe lattices, up to the subtleties discussed above. 

5.2 Existence of frozen phase 

KA models on a Bethe lattice with connectivity k + 1 are defined as in section 
11.31 with < m < k and s = k — m. As usual, it is convenient to arrange the 
lattice as a tree with k branches going up from each site and one going down. 
Before proceeding with a more careful analysis, we first give an argument that on 
Bethe lattices, in contrast to hypercubic lattices, infinite frozen clusters almost 
always exist at sufficiently high densities, p. 

Rough analogs of the infinite networks of frozen fully-occupied slabs that 
can occur on hypercubic lattices with s < d, are fully occupied infinite subtrees 
with branching ratio m = k — s and hence coordination number m + 1: none of 
the particles on such a fully occupied subtree can move. If a Bethe lattice with 
coordination number four is arranged with the four bonds coming out of each 
vertex forming a cross, then all the vertices of a fully occupied three-sub tree 
will appear to have a line going straight through them, and another half line 
that cither terminates at the vertex in a T-junction, or also continues through it. 
Such trees are thus somewhat similar to a network of bars on a square lattice. In 
another sense, however, because they have no loops, they are more similar to a 
single two-wide infinitely-long solid bar on a square lattice. Because, for m < k, 
the number of potential m-sub trees is exponentially large — in contrast to the 
number of infinitely long straight slabs on hypercubic lattices — it is natural 
to expect that a fully occupied frozen subtree will almost always exist at high 
densities. This can be shown by a simple generalization of the analysis for 
conventional site percolation on Bethe lattices [48,55]. 

Consider a site, call it i, and the probability, Q' , that it is both occupied and 
is the root of a fully occupied m— subtree on the portion of the Bethe lattice 
above it. As the existence of such a subtree requires the existence of similar 
subtrees rooted on at least m of the sites above site i, and these are independent 
events, each with some probability Q, we can write a simple recursion relation: 



This has the trivial fixed point solution Q = for any density, but, above a 
critical density, pr, it has an additional non-trivial fixed point solution, Q*(p) > 
0. For the simplest interesting case, k = 3, s = 1 — loosely approximating the 
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coordination-four square lattice — , one can easily solve the fixed point equation 
to find pt = §, a discontinuous jump at p = pr to a non-zero value: Q*(pr) = § 
and above px a square-root cusp increase in Q* (p) which eventually saturates at 
Q*(l) = 1. Other cases with m > 2 behave qualitatively similarly, in particular 
also exhibiting a discontinuous onset of the fraction of sites in the frozen sub- 
tree — simply related to Q* — at a non-trivial critical point at p = pT(k,m). 
For models with m = 1, i.e. s = k — 1, the existence of an m-subtree is the 
same as the existence of an infinite fully occupied cluster and the behavior of Q 
is thus identical to that of conventional percolation: there is a critical point at 
Pt = \ with Q* rising linearly from zero as p increases further. 

What is the significance of these results for KA models? As do other perco- 
lation arguments, these provide useful bounds: the existence of an infinite fully 
occupied m-tree on a /c-Bethe lattice, i.e. Q* > 0, is a sufficient condition for 
the existence of an infinite frozen cluster in the corresponding fc, m KA model. 
Thus we can immediately conclude that a frozen phase is possible in these Bethe 
lattice models at least for p > pt- We do not, however, yet have a necessary 
condition for frozen clusters: as can be seen, these can have isolated vacancies 
within them at least at sites with more than m + 1 frozen subtrees nearby. 



5.3 Dynamical transition 

We now analyze the behavior with the actual KA constraints by a generalization 
of the downwards iterative procedure used above. 

Consider a particular site and define the following events: 

(i) The site is occupied by a particle which can never move up as long as the 

site below it is occupied; denote by Y the probability of this event. 

(ii) The site is frozen: occupied by a particle which can never move up even 

if the site below it is empty; denote by F the probability of this event, 
which is a subset of event (i). 

(iii) The site is empty but blocked in such a way that particles on lower sites 
can never move up to the site; denote by B the probability of this event. 

By using the tree-like structure and its symmetry under exchange of branches, 
one can write iterative equations for the (primed) probabilities for the sites in 
one layer, in terms of the (unprimed) probabilities for the sites in the layer 
immediately above: 
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These can be alternately written for Y in one layer as a function of the three 
Ya in the three next higher layers. 
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We are particularly interested in stable fixed points of the iterative equations 
which correspond to the desired behavior in the interior of the system; the fixed 
point equations can be reduced to independent polynomial equations, with p 
dependent coefficients, for any one of the Y, B, or F. For ease of notation, 
except when it might be confusing, we will denote fixed point probabilities 
simply by Y, F, B and Q. 

We have earlier shown that there exists a frozen phase at high densities: 
in particular for p > px the probability Q is non-zero. Since when the event 
corresponding to Q occurs, the particle is definitely frozen, we know that F > Q; 
F non-zero then implies from 1)59(1 that both Y and B are also non-zero at high 
densities. 

To show that all the restricted events have zero probability at low densities, 
it is sufficient to note that Y', F' and B' all depend at most linearly on Y, F, 
and B (as the other powers in the polynomials are greater than unity) , and two 
of the coefficients are p. Thus at sufficiently low densities, Y must decrease on 
iteration by a factor that is no larger than a multiple of p. Thus the only fixed 
point at low densities is Y = B = F = 0. 

We have thereby shown that a dynamical phase transition must exist for 
Bethe lattice KA models with < s < k. We now show that, except for 
s = k — 1, this transition is discontinuous. 

5.3.1 Discontinuity at transition for s < k — 1 

For KA models with s < k— 1 , i.e.: m > 2, one can show straightforwardly that 
the transition is discontinuous. Focus on one site, call it x. A particle on x can 
move up even if the site below it is occupied, if at most one of the k sites above 
x and the k 2 sites one level further up are occupied. In such a situation, the 
marginal event whose probability is Y does not occur. For the empty sites in 
the two levels above x, the marginal occupied event does not occur either and 
the probability of such sites is thus at least 1 — Y. With all probabilities having 
their fixed point values, these observations imply that 

Y < 1 - (1 - Y) k+k2 -(k + k 2 )Y(l - r) fc + fc2 -! . (60) 

But this inequality cannot be obtained if Y is positive but small, since for 
< Y <C f it would become essentially 1 < (k + k 2 — l)(k + k 2 )Y, which cannot 
hold for small Y. Therefore, the fixed point value of Y cannot go continuously 
to zero: it must jump to a non-zero value at the critical density, p c . 

Note that this argument does not apply for the case s = k — 1. This is 
because, even if all but one of the sites immediately above x are empty, a 
particle on x can still not move up if the site below it is occupied (the inequality 
(I60|l in this case does not have the second term on the right hand side; this 
invalidates the rest of the argument). 

5.3.2 Continuity of transition for s=k-l 

Now consider the most highly constrained non-trivial case: s = k— 1. In contrast 
to the less constrained models, for a particle to be frozen in this case, it need 
only belong to an infinite cluster of particles with two or more neighbors, i.e. to 
a percolating cluster in the conventional sense of nearest neighbor percolation. 
This means that for m = 1, the critical density for the KA model should be 
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the same as for percolation, p c = 1/fc [55] , and the transition will be continuous 
with Y growing linearly just above p c . By solving equations (|59|l for the simplest 
such cases, k = 2 and k = 3, this can be checked directly. More generally, for 
any k > 2 with s = k — 1, to quadratic order when it is small, V depends only 
on Y and can be seen to attain a fixed point value Y rss -j^j (p — p c ) just above 
the transition at p c = 1/fc; this non-zero Y induces non-zero, but quadratically 
small, values for both B and F. 

5.3.3 Connection to bootstrap percolation 

As for hypercubic lattices, bootstrap percolation yields bounds for the behavior 
on Bet he lattices. Before analyzing in detail the character of the dynamical tran- 
sition, we illustrate this connection which can elucidate the physical mechanism 
underlying the dynamical transition of KA models. 

The relevant bootstrap percolation process is similar to that given earlier 
for hypercubic lattices (and introduced by Kob and Andersen [11]): From a 
random initial configuration with density p, remove all the particles that can 
move according to KA rules and iterate the process until no more particles can 
be removed. If some particles survive at the end of this procedure, they must be 
frozen in their initial configuration under the KA dynamics. Define p B as the 
probability that at the end of the procedure an infinite particle cluster remains. 
By the above observation, if p B > the KA model at density p is not ergodic. 
Let p B be the critical density, if any, at which a bootstrap percolation transition 
takes place, namely p B — for p < p B and p B > for p > p B . If there exists 
such a bootstrap transition, a dynamical transition for the corresponding KA 
model must also occur and the corresponding critical density satisfies p c < pb- 
Moreover, a reasonable Ansatz, is that the two transitions and therefore the 
corresponding critical densities, in fact, coincide. In the rest of this section we 
present an argument supporting this conjecture. 

Let a; be a site on which a particle is frozen with respect to KA dynamics, i.e. 
a particle that cannot move on any time scale after the thermodynamic limit 
has been taken. Than x must have either, more than m = z — s — 1 neighbors 
that are occupied and frozen, or, if to or fewer of the neighbors are frozen 
(and thus the particle on x could potentially move), all its unfrozen neighboring 
sites, either empty or occupied, must themselves have more than to occupied 
frozen neighbors. If this were not the case, one could move all the unfrozen 
neighbors out of the way, and then move the putatively frozen particle from 
site x. Moving the neighbors at the same time is possible in such a situation 
because the constraints on them are independent since loops that could induce 
such correlations, do not exist on Bethe lattices,. In other words, moving one 
of the neighboring particles away from a site should not affect the ability of its 
other neighbors to move and an unfrozen neighboring particle can itself move 
provided an appropriate set of particles that are further up its branch of the 
tree have moved out of the way. One can repeat the above argument to show 
that the neighbors of x must each have at least to — 1 = 2 — s — 2 neighbors 
along their own branches occupied by frozen particles, or else all their unblocked 
neighbors must themselves have more than z — s — 1 frozen neighbors, and so 
on and so forth. Therefore, the probability P that a given site is occupied by a 
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frozen particle can be expressed in terms of Y and B as 



(fc+l k+l \ 

j = fc-s+l J=a+1 / 

As a consequence in the ergodic phase Y = B = and P = 0, whereas in the 
non ergodic phase Y and B are strictly positive and thus so is P 13 . 

This result leads to a physical interpretation of the dynamical transition 
as a "jamming" transition: at p c , an infinite set of blocked particles suddenly 
appears and concomitantly a breaking of ergodicity occurs. 

We have checked the relation between the bootstrap percolation transition 
and the dynamical transition by performing numerically the bootstrap proce- 
dure for the case k = 3,s = 1, and verifying that the density above which a 
cluster of particles remains at the end of the process is compatible with the 
critical density for the dynamical transition. 



5.3.4 Quantitative results for k = 3, s = 1 and k = 5, s = 2 models 

We now consider in detail the cases k = 3, s = 1 and k = 5, s = 3, which mimic, 
respectively, the non-trivial square lattice model, and the three-dimensional 
cubic-lattice model originally introduced by Kob and Andersen. 
In the case k = 2 s = 1, equations (|59|) give 



Y' = p(Y 3 + 3F 2 (1 -Y) + £ 3 + 3B 2 Y) 
B> = (l-p)F 3 

F' = p{Y + Bf (62) 

from which one can immediately see that, at a fixed point with Y > also 
B > and F > 0. Therefore one can study the transition by analyzing the 
fixed point equations for 

G = Y + B (63) 
which can be written in polynomial form: 



R(G) = G[-l + 3pG-2pG 2 + vp 3 G s -6vp i {G 9 -G w ) + 3v 2 p 7 {G 17 -G 18 )] = . 

(64) 

By direct analysis one can see that there is a fixed point with G = which 
is stable for any density, while above a critical density, p c , a second physical 
solution appears with positive G; both are stable in this regime. Separating 
these stable fixed points is an unstable one at an intermediate value of G. As 
the density decreases to p = p c , this unstable fixed point annihilates with the 
non-zero stable fixed point at a saddle-node bifurcation signaled by 

2S|M = (65) 



13 We expect that our derivation of P is correct under the hypothesis that to pick away 
a particle during the bootstrap procedure, it is not necessary to pick away before so doing, 
an infinite (i.e. diverging with the system size) number of other particles. This should be 
roughly equivalent to the hypothesis that an appropriately defined correlation length is finite, 
a reasonable expectation except at the bootstrap transition point. 
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Figure 10: The probability, G (defined in the text) as a function of density for 
k = 3, s = 1 Bethe lattice as obtained from the stable solution of equations (|64|) . 
In the inset log(G — G c ) is shown as a function of log(p — p c ) in the vicinity of 
p c showing the square root dependence. 

at the fixed point. 

The values of the critical density p c and critical G c = G(p c ) are: 

p c ~ 0.888 G c ^ 0.758 (66) 

The critical density is strictly less than, although only very slightly so. the 
critical density pT = § , of the much simpler two-branching-subtree percolation 
problem discussed in section i|5.2p . 

At p = p c the value of G jumps discontinuously from the low density value 
G = to G c and then increases with a square root cusp: G « G c + C^p — p~ c j p c , 
see figureEH Similar behaviour obtains for the probabilities P of a site occupied 
by a frozen particle, as can be obtained from equations H61|) and (|62|l . 

For the case k = 5, s = 2, by solving equations l|59|) . a similar transition is 
found at a slightly higher critical density p c ~ 0.915 

5.3.5 Diverging time and length scales at the dynamical transition 

Although the dynamical transition in these KA models is discontinuous, there 
is precursor behavior more characteristic of critical transitions. In this respect, 
the situation is somewhat analogous to conventional bond percolation on one 
dimensional lattices with connection probability between any two points pro- 
portional to the inverse square of the distance between them [56] . 

Indeed from the discussion in the previous subsection, an interpretation of 
the mixed nature of the dynamical transition of Bcthe lattice KA models is 
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suggested: the dense cluster of frozen particles that arises immediately at p = p c 
contains a finite fraction of the particles and it is, at the same time, fragile in 
a way similar to incipient infinite clusters at conventional critical percolation. 
This fragility implies that by removing a tiny, but carefully chosen, fraction of 
the particles, the frozen cluster can be "melted" and will disappear. Associated 
with this fragility is a correlation length that diverges as p — > p c from above. 
The form of the divergence, ~ (p — p c )~~^ is controlled by the "rate" - in tree 
levels - of approach of the iterated probabilities to the fixed point. Physically, 
this correlation length is the (chemical) distance beyond which it is unlikely to 
have to move vacancies in order for a given particle to be moved. 

More significant than the long length scales that appear near the critical 
density, are the corresponding dynamical effects. In particular, we expect a 
diverging time scale as p c is approached from below. Unfortunately, our analytic 
methods, while they take into account the dynamical constraints, cannot directly 
be used to compute the dynamical correlations of primary interest. At this point, 
we resort to numerical simulations to study these and other properties of the 
Bcthe lattice models. 

Before turning to the numerics, we note that the behavior on Cayley trees, 
because of their boundaries which account for a positive fraction of the volume, 
might exhibit different phase diagrams depending on the choice of boundary 
conditions. 



5.4 Dynamics and numerical simulations 

In order to understand in more detail the character of the dynamical transition, 
we have performed numerical simulations of KA dynamics for the k = 3 s = 1, 
model on a random c-regular graph of N = 10 4 sites with coordination number 
k + 1 = 4. 

In particular we have computed the local density-density dynamic auto- 
correlation function C(t) and the corresponding dynamical susceptibility 
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where ( ) denotes an average over the Bernoulli product measure at density 
p for the choice of the initial configuration at over the Montecarlo dynamics. 
The dynamical susceptibility Xi(t) (and its generalization) has been introduced 
in [57] and is currently used [7, 24] to detect a growing dynamical correlation 
length in simulations of glass- forming liquids. 

The results for C (t) , plotted in figure ^2 show that the equilibration time 
grows and seems to diverge at the critical density found from the quasistatic 
analytical calculations (see equation JSBJ). Indeed, for p < p c — A, C(t) tends to 
zero at long times but with a characteristic decay time that increases at higher 
densities. Here A, the precision obtained for p c , is 0.0075. For p < p c + A, 
C(t) no longer appears to decay to zero but seems to approach a non-zero 
plateau whose magnitude it is natural to call an Edwards- Anderson parameter 
by analogy with spin glasses. . Since at the critical density there should be 
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Figure 11: Dynamic autocorrelation function C(t) as a function of time for den- 
sities 0.85, 0.86, 0.87, 0.875, 0.9 (from down to up) on a 10000-site approximation 
to a Bethe lattice with coordination number k + 1 = 4 and single-vacancy as- 
sisted dynamics. The dynamical transition is clearly evident between the den- 
sities 0.875 and 0.9. The straight line is the value of the Edwards-Anderson 
order-parameter obtained by the approximation discussed in the text. 

NP(p c ) frozen particles, the value of the Edwards Anderson parameter qEA = 
lim^oo limpf^oo C(t) has two contributions coming from sites that are occupied 
respectively by a frozen particle or by a vacancy that can never move. These 
can be computed exactly from the quasi-static computation developed in the 
previous section. However there is a third contribution that corresponds to 



where the index i = 1, . . . , N' runs only over the sites whose occupation num- 
ber changes during the dynamics (i.e. they are not occupied by vacancies or 
frozen particles). This contribution to qEA is more cumbersome to compute. 
In this paper we just give a rough approximation which consists of assuming 
that the long-time limit of (ni(t)rii(0)) on these sites — which consist of many 
small disconnected clusters — is equal to the density of particles within the full 
set of N' dynamic sites. This can be computed from the quasi-static analysis 
developed in the previous section. Comparison with numerics shows that this 
approximation works rather well (see e.g. Fig. II ip . 

A dynamical quantity that can be computed exactly is the long-time limit 
of the persistence function: the probability that the occupation variable of a 
given sites has not changed from its initial value for an interval of time t. The 
long-time plateau of this, by definition, is equal to the fraction of sites that are 
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Figure 12: Dynamical susceptibility Xi{t) as a function of time for densities 
0.85,0.86,0.87,0.875 (from down to up) on a 10000-sitc approximation to a 
Bethe lattice with coordination number k + 1 = 4 and single- vacancy assisted 
dynamics. 

occupied by either a frozen particle or by a vacancy that can never move. 

Plots of the dynamical susceptibility, X4(t), found in the simulations are 
shown in figure IT2"1 for various densities: these appear consistent with a diver- 
gence at the critical density. More precisely both the height of the peak and the 
characteristic time scales at which this occurs, appear to diverge at p c . 

5.4.1 Conflgurational entropy 

If phase space breaks into a sufficiently large number of different ergodic com- 
ponents, some of the equilibrium entropy density will be associated with the 
logarithm of the number, Af e , of distinct components: 

S c {p) = hxM e /V (67) 

where V is the total number of sites. This is referred to in the glass literature 
as the conflgurational entropy 14 . 

For KA models on Bethe lattices one can find a lower bound for the conflgu- 
rational entropy by observing that all the configurations belonging to the same 
ergodic component must have the same cluster of frozen particles, therefore Af e 
is at least as large as the number of different sets of frozen particles, Aff, which 

14 This terminology is somewhat misleading since in general statistical mechanics "conflg- 
urational entropy" refers to the sum over all configurations (after having integrated out the 
momenta). For glasses, its meaning is the contribution to the entropy density coming from 
the large number of different ergodic components. The other contribution to the entropy is 
related to the number of configurations within the given ergodic component. 
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Figure 13: Lower bound, S b (p), on the configurational entropy — logarithm of 
the number of frozen configurations per site — as a function of density for the 
k — 3 s = 1 Bcthc lattice is indicated by the dotted line; this bound is obtained 
from equation l|71|l and numerical solution of equations (|61fl . (|62|) . Below the 
critical density p c ~ 0.888, the lower bound jumps to zero. The straight line is 
the equilibrium entropy. 

has an entropy density, 

S i (p >PF ) = llnjV>. (68) 

Furthermore, the entropy of the non-frozen sites, S t (p 1 pF), can be bounded 
above by the number of ways of putting V p — VpF particles on V — V ps sites, 
where pf is the density of frozen particles at total density p: 

S 1 (p,pb) < {\-pb)H^-Pb)-(p-Pb)Hp-Pb)-{1-p)H^-p) ■ (69) 
Since the total entropy density, 

S(p) = -p\np- (1 - p)\n(l - p) = {Sf(p,p~ F ) + S c (p)) , (70) 

and we know that p~F = P that was computed earlier, see eq. (|fj 1|> . we obtain a 
lower bound for the configurational entropy 

S c (p) > S f (p,P) > -p\np-(l-P)\n(l-P) + (p-P)\n(p-P) . (71) 

In figure 1131 this lower bound is plotted as a function of density for the case 
k = 3, s = 1. Note that, among other factors, we have ignored in this estimate 
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Figure 14: A branch of a Bethe lattice with triangular loops. 

the existence of empty but blocked sites: the effect of these will be to lower the 
entropy of the unfrozen particles. 

In conclusion, our results from lower bounds and analytical and numerical 
analysis of the dynamical transition support the expectation that the config- 
urational entropy jumps form zero to a non-zero finite value precisely at p c . 
For p > p c , the configurational entropy decreases with the decreasing vacancy 
density and overall entropy; it vanishes at p = 1. 

5.5 Bethe lattices with loops 

In order to investigate whether the absence of loops is the essential reason for the 
existence of a phase transition on Bethe lattices, we briefly consider structures 
that are tree-like on large scales, but do have loops. 

Specifically, consider the rooted tree composed of triangles with one vertex 
pointing downwards and the other two vertices each being the bottom vertex of 
another triangle and hence the root of a tree. Then define a Bethe lattice with 
triangular loops as the graph obtained by taking two such lattices and merging 
the free vertices of the two roots (see figure I14|l . The coordination number of 
such a tree is z = 4 and the branching ratio k = 2. This is an example of a 
cactus, or Husimi, tree [58] with primary clusters of three vertices and three 
branches departing from each cluster. 

Consider the KA model with s = 1 on the triangular tree defined above. It is 
useful to focus on the set £{ of configurations of the subtree rooted at triangle i, 
such that the bottom vertex is occupied and a vacancy cannot be brought down 
to it without using vacancies below it on the tree. Let P be the probability 
of such an event. If a configuration, r\ has all three sites of triangle i occupied 
and all of the subtrees rooted on triangles two levels up from i, i 1 , . . . , i 4 , have 
configurations £ then the configuration of the subtree rooted at triangle i 
belongs to £. L . The same is true for any configuration, 77, that has the three sites 
of i all occupied, and r) € £y for three of the four subtrees rooted at the {V}- 
Therefore 

P> f(P) = P 3 {P 4 + 4:P 3 {1- P)) . (72) 

Since P = 1 is a fixed point of P' = f(P) when p = 1 , at sufficiently large 
density there exists a P > such that f(P) = P. If one iterates down the 
tree for actual Ps, then the inequality (|72|l implies that, starting from an initial 
condition P = P + e > P, all the subsequent P's, Pi,P%, . ■ . , will be larger 
than P. Indeed, since the f(x) is monotonically increasing in < x < 1, 
Pi > f(P + e ) > f{P) = P an d the same is true at any subsequent step. 
Therefore, the full fixed point must have P* > P > 0. Thus for sufficiently 
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large p, P*{p) will be non zero and the system in a frozen phase. It immediately 
follows that the density of frozen particles, P is non-zero. 

A generalization of the argument in section 15.21 allows one to prove that 
there exists a transition at a non-zero density from a regime where P = to a 
region where P is positive. Let £\ be the complement of defined above. A 
configuration is definitely in £ j if the four branches rooted at the four triangles 
two levels up from i are in . The same is true if just three of them are in 
£ p and furthermore one can have or bring vacancies in the highest sites of the 
triangle which is just above the site which is in the same triangle of the site 
that cannot have the vacancy. This sufficient condition for £} yields and upper 
bound on P: 

P< 1-(1-P) 4 -4P(I-P) 4 (73) 

which implies — as in section l5.3.1l — that P has a jump from zero to a non-zero 
value at the critical density. 

The arguments given above can be extended to KA models with any con- 
straint parameter, s, on generalized Bethe lattices with n-polygonal loops for 
any n, i.e. trees composed of polygons of n vertices, each which is also a vertex 
of another polygon. To prove the existence of the transition, it is sufficient to 
consider the set of configurations £j such that a vacancy is not and cannot be 
brought down to the bottom vertex i without using vacancies below. By the 
same argument as above, we can bound the probability P of such an event by 

P > p n (V"- 1 ) 2 + (n - i)2p(»-i) 2 -i(! _ P f\ ( 74 ) 

which implies again that P > for sufficiently high density. Generalizing the 
procedure used for the triangular tree one obtains an bound in the opposite 
direction: 

P < 1 - (1 - P)^ 1 - 1 ) 2 - ( n - ifp(^) 2 -\i _ p) (75) 

which implies, as above, that the transition is discontinuous. 

Another class of Bethe-like lattices with loops is also instructive to analyze. 
Consider, as an example, starting with normal Bcthe lattice with connectivity 
four and replacing each site by a square lattice A G Z 2 of linear size L connected 
to to other squares along its edges (see figure IT5|l . If we arrange the lattice as a 
tree with one branch pointing down from each vertex, the "upwards" branches 
from a square consist of an upwards extension of the lower branch and side 
branches emerging perpendicularly from the other two sides. Locally, this looks 
like a square lattice except at the corners: the corner sites still have four neigh- 
bors — albeit on three other squares — but have eight second neighbors that are 
each on separate "sheets" in contrast to the four second neighbors on a regular 
square lattice each of which can be reached in two distinct ways. 

Drawing the tree of squares as described above immediately implies that for 
the s = 1 KA model on which we will focus, there exist infinite networks of 
frozen particles. As for the normal square lattice, these consist of two-wide fully 
occupied bars, each either extending to a faraway boundary, or beginning and 
terminating at T-junctions with other such bars. But such a network, because 
of the large scale topology of the Bethe lattice, is very different from its real- 
square- lattice counterpart. It is still very fragile but there are so many potential 
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Figure 15: Bethe lattices with loops (c) obtained by substituting sites with 
squares of size L (b) from a normal Bethe lattices with connectivity four (a). 
Here L = 3. 



T-junctions, which are independent of each other if the bars emerging from them 
go onto different branches, that such a frozen network exists almost surely at 
sufficiently high densities. The reason for this is simple: consider a vertical 
bar coming into a square, i, from below. There are at least two possibilities: 
the bar can continue through i and into the square above it, or it can end 
at a T-junction with a bar that extends horizontally into both of the other two 
neighboring squares of i. Since the probabilities, T, that the extensions into each 
of the three neighboring squares exist and arc themselves part of infinite frozen 
networks on their own sub-trees — i.e., excluding square i — are independent 
of each other, one can bound from below the probability that square i has a bar 
rising vertically from its lower boundary and either extending through the top 
boundary or ending at a T-junction with a bar that extends into the two side 
boundaries, in terms of T and p. As in cases we have already analyzed, this 
can be used to show that there must be an infinite frozen network at sufficiently 
high densities. One can likewise show that the transition to the frozen phase 
will be discontinuous. 

To show that the transition on this Bethe lattice of squares is discontinuous, 
it is useful to focus on the framing by vacancies of individual squares, in the 
sense used for the conventional square lattice. Consider the event that, without 
taking advantage of possible vacancies in the square below /, one cannot frame 
the configuration in /, i.e. reach a configuration with all the boundary sites 
of / empty. Call X the probability of such an event. If, without using the 
square i, the three squares left, right and above I can all framed, or if either 
the above and left or the above and right squares can be framed, or if the left 
and right squares can be framed and in the square / the first row is empty and 
all the others are occupied then the square / can definitely be framed. These 
possibilities are sufficient, but not necessary. Therefore 

X < 1 - (1 - X) 3 - [2 + p L2 ~ L (l - p) L ]X(l - Xf (76) 

which is incompatible with a continuous transition by an argument similar to 
that used for other Bethe lattices. 

In general, the existence of finite loops on these decorated Bethe lattices, 
makes the motion of particles less restricted that on a regular Bethe lattice. 
Indeed, there will be configurations in which two or more vacancies can be 
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ejected from one branch up into another, and by making use of these, trigger 
the potential ejection of a larger number of vacancies down from the second 
branch. This process can depend on configurations arbitrarily far up in the 
two branches. Such ejected vacancies could in principle extend their influence 
arbitrarily further down the tree and render all particles mobile, even when the 
density is above the critical value for the loopless Bethe lattice. But the above 
results show that for low enough vacancy concentration this will not occur: the 
system will be in a partially frozen state like that of the simple Bethe lattices. 
The above picture does, however, suggest a route to estimating the critical 
density as a function of L. 

The trees-of-squares we have been studying interpolate between simple Bethe 
lattices, for L = 1, and a full square lattice for L = oo. Therefore on the trees- 
of-squares, in the limit L — > oo the critical density, p c (L) should go to one. 
But how does it approach one? From scaling arguments, one would guess that 
the framing probabilities on the tree will be functions of the ratio of L to one 
of the characteristic lengths of the infinite square system: i.e. functions of 
either L/£(p), or Naively, one might expect that the behavior would be 

controlled by the large crossover length, which determines the spatial density, 
1/S(p) 2 , of mobile cores and thus the probability that squares of size L are 
frameable. However, this ignores the effects of possible help from vacancies 
moved from other squares and the correct behavior can be seen as follows. 
Consider one of the L by L squares that can be framed by vacancies. The 
probability that this framing can be continued all the way through a given one 
of its neighboring squares to make an L x 2L framed rectangle, is of order 
[a(l — p)L] L with an order-one coefficient a. This will be substantial as long 
as £ (p) > L — in spite of the fact that the probability the original square was 
frameable is very small unless L ~ The crucial observation that leads 

to a scaling with £, is a consequence of the tree structure: once the neighbor 
square of i has been framed together with i, as above, each of the three other 
neighbors of this neighbor can potentially join the frame as there already exist 
complete rows of vacancies along one of their sides from which the framing can 
be extended. Since any of these squares that joins the frame gives rises to 
another three potential squares that might join as well, if the original square 
frame can be extended out in several directions to make a small frameable tree, 
it become likely that it can be extended much further, indeed, to infinity. Thus 
the crucial condition for extending the framing to infinity is that there be a 
substantial probability of a frame around one square being extendable across 
one of its neighboring squares — even when the probability that a single square 
is frameable on its own is beery small. 

We conclude that the critical density, p c (L) for the s = 1 KA model on a 
tree of L x L squares, will be given by 

^(Pc) = bL (77) 

with some order-one coefficient b. As £ ~ ln[l/(l — p)]/(l — p) for high densities, 
this yields 

BhxL 

Pc (L) » 1 - — — (78) 

with some coefficient B. Note that this form is the same as that suggested by 
the argument given above for the appearance of an infinite network of bars of 
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frozen particles, since the probability of a two-wide fully occupied bar crossing 
a square is of order Lp 2L which is of order one at p c (L). 



6 Fredrickson Andersen model 

The Fredrickson Andersen (FA) model is a spin model with kinetically con- 
strained dynamics that is similar to the KA model but without a conservation 
law: roughly speaking KA models are FA models with dynamics that con- 
serves the number of particles. Some of the results discussed previously for KA 
models have been previously obtained for FA models [9,28,60]: for example 
irrcducibility on some hypcrcubic lattices for the non-trivial range of the con- 
straint parameter . However, there are important results, like the absence of 
a crgodic/non-ergodic transition, that have not been proven previously, as well 
others, such as the dependence of the relaxation time scale on density, that can 
be made sharper. In the following we briefly discuss how these new results on 
FA models can be obtained using the same techniques that we have developed 
for KA models. 



6.1 Definition of the model 

Let A be an hypcrcubic d-dimensional lattice and / an integer parameter in the 
range 0, ..2d. The FA model [28] is a facilitated Ising model with "occupation" 
variables r\ x £ {0, 1} — configuration space 17a = {0, 1}' A '. The FA model is 
conventionally defined in terms of spin variables; we use occupation variables 
to make comparisons with KA models more natural.] The occupation variables 
can "flip" with rate 

I \ - / if Ux ^ - f (7Q\ 

Cx{T,) — \ min(l,exp-^2) otherwise l ' yj 

where 

n x (v)-= E C 1 -^) > ( 8 °) 
2/eA 

\x-y\=l 

i.e. n x is the number of empty neighbors of x and the change in energy AH = 
H(rf) — H(rj) comes from the non interacting Hamiltonian 

X 

with h a positive coefficient that determines the fraction of occupied sites in 
equilibrium. These rates satisfy detailed balance with respect to the trivial 
equilibrium product measure 

with p = 1/T: this corresponds to Bernoulli measure with density 

p = exp(/3h)/(l + cxp(/3/i)) . (83) 
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In contrast to KA models, Ylx£A^ x is not conserved in FA models: a "flip" at 
a single site corresponds to the birth or death of a particle. 

In the unrestricted case / = 0, the FA dynamics is simply Glauber dynamics 
for uncoupled Ising spins in a magnetic field: occupied sites corresponding to 
up spins and vacant sites to down spins. But for positive /, a site can change 
only if at least / of its nearest neighbors arc empty (i.e. their occupation 
variables equal zero), we refer to empty sites as facilitating and to / as the 
facilitation parameter. For / > d, fully occupied hypercubes will be frozen. 
We thus restrict consideration to 1 < / < d. Since energy favors occupied 
sites, in the low temperature regime most potential flips will not be allowed. 
Therefore we expect dynamics to be slow in this regime, as for the high density 
regime of KA models. On finite lattices, as for KA models, there is not a unique 
invariant measure: there exist blocked configurations and the FA dynamics are 
not ergodic. But irreducibility in the thermodynamic limit for / < d has been 
proved for two and three-dimensional hypercubic lattice FA models, [49,60,61] 
(see [9] for a review). 

There is considerable similarity between KA and FA models: in the former, 
a particle cannot move if it has more than m occupied neighbors; in the latter 
it cannot flip to vacant if it has more than 2d — / occupied neighbors. Thus 
/ plays a similar role to 2d — m. or, equivalcntly, to s + 1, with s the KA 
model parameter that restricts a vacancy from moving if it has fewer than s 
neighboring vacancies. In addition to their conservation law, the KA models 
differ from FA models because of their restrictions on the neighbors of two sites 
(the initial and final sites for the jumping particle) for each possible move. 

6.2 Irreducibility and ergodicity 

The presence or absence of frozen configurations in FA models is simply related 
to bootstrap percolation. Consider a configuration of occupation variables and 
perform a boostrap procedure by removing at each step all particles that have 
at least / empty neighbors. If, at the end of this procedure no particles remain, 
the initial configuration belongs to the irreducible component that contains the 
totally empty configuration (called the high temperature partition by Fredrick- 
son and Andersen). This is exactly conventional bootstrap percolation. But all 
the "moves" made in this bootstrap procedure are allowed in the correspond- 
ing FA model. Therefore, using the results for the absence of an unpercolatcd 
phase for bootstrap percolation [49,61], one can immediately conclude that in 
the thermodynamic limit irreducibility holds for FA models at any tempera- 
ture as long as / < d. In contrast to KA models for which bootstrap moves 
would not be allowed, for FA models bootstrap results are sufficient to conclude 
irreducibility in the thermodynamic limit. The result proven for the two and 
three-dimensional cases [60] is trivial to generalize to higher dimensional cases 
using the bootstrap percolation results [62]. Furthermore, the crossover length 
L B (p) that characterizes the bootstrap procedure is not only a lower bound for 
the crossover length, S(T) of FA models, (as was the case for KA models) but 
coincides with it under the change of variables p — > exp(/3ft,)/(l + exp(/3/i)). 

As we explained previously, irreducibility in the thermodynamic limit is not 
a sufficient condition for ergodicity, which is the physically interesting property 
for dynamics. However, ergodicity can be proven as in section \2. 51 using again 
irreducibility and the product form of equilibrium measure. 
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6.2.1 Transition on Bethe lattices 

On Bethc lattices, FA models with positive / displays a dynamical transition 
and aging dynamics [59,63]. This is a natural expectation that follows from 
the bootstrap percolation transition on Bethe lattices: our analysis for the KA 
model — part of which involved removing particles as allowed in FA models 
- leads, in the FA case on a Bethe lattice, to exactly the same bootstrap 
percolation equations as in [48]. We will not discuss further Bethe lattice FA 
models. 

6.3 Dynamics 

We now turn to the dynamics for FA models on hypcrcubic d-dimcnsional lat- 
tices. 

It is well known from previous results and simulations that the behaviour is 
strikingly different for the case / = 1 than for larger / [5,9,28]. With / = 1, a 
single vacancy can facilitate the flip of any of its neighbors, therefore relaxation 
does not require cooperative processes and time scales are proportional to the 
density of vacancies, as confirmed by numerics. This least-restricted FA model 
is loosely analogous to the s — 1 KA models normal hard core lattice gas (which 
would correspond to KA with s = 0) for which a single vacancy can freely move 
in an otherwise totally filled lattice: in the FA case, a vacancy can "move" by 
flipping a neighboring occupied site, and then flipping to occupied at the original 
site. Nevertheless, this analogy cannot be take too far: interesting — although 
not cooperative — behavior does exist for / = 1 FA models. Indeed it has been 
shown in [19] that in dimension less than four vacancies cannot be considered 
as independently diffusing defects (in contrast to the trivial KA model in which 
single vacancies can move). The interaction between vacancies, in particular 
their creation and annihilation, lead to a reaction-diffusion process that belongs 
to the same universality class as directed percolation in three dimensions. 

For more restricted FA models with the facilitating parameter / > 1 neither 
a single facilitating vacancy, nor a finite cluster of them, can enable all the other 
sites to flip if these are initially all occupied: as for KA models, a finite vacancy 
cluster can never breakthrough a slab that is infinite in 2d— f + 1 directions and 
of width two in the other / — 1 directions. Hence, cooperative processes must 
be involved in the low temperature relaxation of these FA models. The nature 
of the cooperative mechanism that causes equilibration has been conjectured by 
Reiter [60]: it is very similar to the one we discussed for the KA model in [37] 
with the correspondence / — ► s + 1. 

Consider for simplicity the two-dimensional case with / = 2 at temperature 
T and focus on regions of an infinite lattice that are frameable — - in the KA 
sense with s = 1 (see section l2~T)l — - out to linear size £(p), with p(T) defined in 
(|5"2l and £ as in (j5J). These we refer to as mobile cores. [37] It is easy to check 
that using moves allowed by FA rules every occupation variable inside the core 
can be changed to zero. On a typical line segment of length £(p) outside the 
special core regions, there is likely to be at least one vacancy. The existence 
of such vacancies around most places a core could be, ensure that a single core 
can diffuse through typical regions of the system — and eventually everywhere. 
Reiter conjectured that the diffusion of these macro-defects is the dominant 
mechanism for relaxation in the low temperature regime of FA models. From 
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knowledge of the temperature dependence of £ and of the spatial density of 
mobile cores, po, — both given by those of the corresponding KA model at the 
same vacancy density, p, see scction lfTTl — together with the characteristic time 
To, within a core, one can make a prediction for the temperature dependence 
of the overall relaxation time. Neglecting, as for KA models, the interactions, 
annihilation and creation of the macro-defects, the bulk relaxation time is [37]: 

. (84) 

PD 

This is the time after which typical sites will be changed because a macro-defect 
has passed by. 

Reiter assumed that the typical timcscalc td for the motion of a macro-defect 
of size / grows as expKl with K a positive constant. However, we have found 
— analytically and numerically — that this is not the case for the KA model: 
Td increases more slowly with I as expi^v^ [37] (see also section0J. This result 
should carry over to the FA case since it is easier for a macro-defect to move 
with FA dynamics than with KA dynamics, yet the dynamics within a core is 
dominated by entropic barriers [37] which will be the same for both models. As 
a consequence of the sub-exponential slowing down within a core, the leading 
contribution to the structural relaxation timescale r is given by the density of 
defects and not by their timescale to move. The same argument should hold 
for any d, 1 < / < d giving r cx p D ~ l/E d j_i(p(T)) d , where Ed,f-i{p) is the 
crossover length of the corresponding KA model and p(T) is the equilibrium 
density of vacancies at temperature T, see eq. (|83|) . For the case d = 2 / = 2, 
we thus have the exact leading low temperature behavior of logr: 

t ~ exp — C °° (85) 
1 - P 

where 2coo ~ 1.1 (see equation l|54f) '). Numerical simulations on a square lattice 
agree well with the above prediction [25] . However one has to be careful because 
of strong finite size effects as discussed for the KA case. And there is another 
complication that might arise in the FA case: the fact that the independent 
defect diffusion assumption fails in the / = 1 case for d < 4 [19], suggests it will 
also for larger /. However if, as for / = 1, the fluctuations in d < 4 just change 
(slightly) the dynamical exponent that relates time and length, then the scaling 
of t with the density for the d — 2 / = 2 case would be the same as in (|85|) but 
perhaps with a (slightly) rcnormalized constant c^. [In the more cooperative 
cases for which there are at least two iterated exponentials, a change in the 
dynamical exponent may not even change the leading asymptotic behavior of 
the relaxation time, but this needs a more careful investigation.] 



7 Conclusions 

We have analyzed in this paper a class of kinetically constrained lattice models 
which had previously been studied primarily via numerical simulations. For the 
KA models on hypercubic lattices, we have found that the behavior in infinite 
systems, in contrast to what was suggested by the simulations, is simple: either 
the system is non-ergodic for all densities because of the existence of finite 
frozen clusters, or it is crgodic at all densities. Nevertheless, in the latter cases 
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dynamics at high densities is intrinsically collective and thus non-trivial. We 
found that there are two length scales that characterize this collective dynamics: 
The smaller one, £, is the linear dimension of a minimal set of vacancies that 
is typically mobile in a system of density p; such regions of size £ d are the 
mobile cores whose motion allows rearrangements of other regions. The longer 
length scale, H, is the typical spacing between mobile cores. Therefore S is 
the crossover length that determines whether finite-size systems will typically 
be almost crgodic, with one irreducible set dominating configuration space, or 
have their configuration space broken up into many large sets. Furthermore, 
the same length controls the high density behaviour of relaxation times, which 
scale as the density of mobile cores, therefore as The manner in which 

these lengths diverge as the density approaches unity depends on the dimension 
and the constraint parameter, s. But, although we have not discussed it thus 
far, there is considerable universality. 



7.1 Universality 

For simple lattices other than hypercubic, the behavior can be guessed by anal- 
ogy with the hypercubic cases. 

If there are no finite frozen clusters, but infinite frozen configurations do 
exist, then the minimum dimension of frozen structures will determine the be- 
havior. If the minimal frozen structures are d — 1 dimensional, then the frames 
of vacancies used to show irreducibility will be one-dimensional and the core 
size will grow with density as, 

ln[l/(l-p)] 

~ (1 _ p) p/( d -i) ( 86 ) 

where p is the number of nearest neighbor vacancies needed on each side to 
enable a frame to grow. The crossover scale grows roughly as S ~ e^. This 
least-restricted collective behavior obtains for s = 1 on hypercubic lattices in 
all dimensions, and also for, e.g., two-dimensional triangular lattices with s = 2 
(s = 1 has mobile pairs of vacancies); and three-dimensional fee lattices (close- 
packed) with s = 3 (s = 1 or 2 models have small mobile clusters of vacancies). 
In all these cases p = 1. 

If the smallest frozen structures are d — 2 dimensional, framing by two- 
dimensional planes is needed: this gives rise to 



e ~ cxp [ (i-py/w ) (87) 

This form cannot occur in two dimensions. In three dimensions, it obtains for 
cubic-lattices with s = 2, as well as for fee lattices, for s = 4 and s = 5, these 
two cases only differing by the numerical coefficient. 

In three dimensions, the above are the only possibilities in the same gen- 
eral class that we have been studying, but in higher dimensions, the iterated 
exponential forms of (|35fl can occur. As we have shown, on Bethe lattices — 
expected to be loosely related to some high dimensional limit — an actual dy- 
namical transition at a non-trivial density, p c , does occur for KA models. 

This immediately suggests a crucial question: Do some finite dimensional 
models exist with this type of dynamical transition? Two candidates in the 
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literature are infinitely thin hard rods on a lattice [64] and a space filling self- 
avoiding chain on a lattice [65]. We leave further exploration for the future. 

But the issue of apparent transitions does already arise in the simple KA 
models. As mentioned at the end of section 0] dramatic slowing down in fi- 
nite dimensions can occur near an almost sharp avoided transition associated 
with a crossover from non-cooperative to cooperative dynamics. In particular, 
some high dimensional lattices — including already some three dimensional lat- 
tices — are approximately tree-like out to some number of neighbors; if the 
dimension increases while the coordination number remains fixed, such tree-like 
local structure suggests that a ghost of the Bcthe lattice transition with the 
same local structure will become more and more apparent. More generally, it 
is likely that in the limit of high dimensions, models that have low-dimensional 
frozen structures — such as s = d — 1 on hypercubic lattices which has frozen 
one-dimensional bars — will, as the density increases, have increasingly rapid 
crossovers that appear more and more like an actual transition. 

In this paper we have studied Bethe lattices with loops that interpolate be- 
tween the finite dimensional lattice and the standard Bcthe lattice. In this case 
there is always a transition although it is pushed to a density that approaches 
one when the finite-dimensional regions become large. From the other direc- 
tion, it would be interesting to analyze a Kac-like interpolation for a fixed finite 
dimensional lattice with longer and longer range connections which converges 
in a certain limit to a Bethe lattice. In this case the transition will presumably 
never occur, but the cross-over will become sharper and sharper near the critical 
density on the Bethe lattice. 

7.2 Comparison with mean-field approaches 

We have found that the generic behavior of KA models on Bethe lattices (excep- 
tions being s = k — 1), has a transition with features similar to both first order 
and critical transitions. Although the density of frozen particles jumps discon- 
tinuously at the critical density, there are precursor effects as the transition is 
approached from the frozen phase and, at least for some quantities, also from 
the ergodic phase. In particular, we have seen that divergence of characteristic 
length and time scales is apparent in the dynamic density-density correlations 
and the related dynamical susceptibility. 

This mixed behavior, as well as several more specific features, is similar 
to that found in certain mean-field-likc theoretical analyses conjectured to be 
relevant for glass transitions. In particular, spin models with quenched random 
interactions that couple collections of p > 2 spins have been studied [66,67]. It 
turns out that the infinite-range versions of such random p-spin models have 
a real thermodynamic transition at a critical temperature, however interesting 
dynamical behavior commences sharply at a higher temperature, To, without 
singularities in static properties. This is interpreted as the onset of ergodicity 
breaking: below Tp the equilibrium measure is fractured in exponentially many 
metastable "states" ("TAP" states which are locally stable solutions of the 
self-consistent mean-field equations); the large number of these give rise to a 
configurational entropy density, S c , that jumps to a non-zero value at To- As the 
temperature is lowered further, S c decreases, vanishing at the thermodynamic 
transition. Due to the presence of quenched random disorder p-spin models 
would seem to have little to do with glasses which must generate their own 
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randomness, but it has been shown that some mean-field approximations to 
non-random models, in particular the mode-coupling approximation, [29] give 
rise to self consistent integral equations with very similar structure to the one 
of this models [67]. In particular, mode coupling equations give rise to an 
analogous ergodicity breaking transition at a finite. We find that the same 
properties are shared by KA models on Bcthc lattices, for which at a finite 
critical density the configuration space is broken into many ergodic components 
and the configurational entropy jumps to a non-zero value and then decreases 
at higher densities (here the thermodynamics is trivial and, hence, Sc equals 
zero only at unit density). In addition, the autocorrelation function, C(t), and 
the associated susceptibility, X4(i), which we computed numerically for Bcthc 
lattices with k = 3 s = 1, behave similarly to those found for the infinite range 
random p-spin models and in the Mode Coupling Theory [57,69,70]. Indeed, 
even the off-equilibrium aging behavior of KCM on Bethe lattices is quite similar 
to the one obtained for mean- field quenched random systems [63] . 

While tempting to conclude that these behaviors are all manifestations of 
the same physical phenomena, one must be careful. In particular the dynamical 
transition in KA models on Bcthc lattices is a reducible/irreducible transition 
whereas for p-spin models the irrcducibility is guaranteed and the breaking of the 
ergodicity is due to the thermodynamic limit. Nevertheless, to explore whether, 
and in what sense, some of the mean-field approximations might capture some 
of the essential physics of more realistic finite dimensional models, the possible 
connections to the KA models are worth exploring further. 

7.3 Experimental issues and prospects 

Kinctically constrained lattice gas models, while they may capture some of 
the important features of glass transitions — or almost transitions — , clearly 
have major flaws. First, in any real system, rates of local processes cannot 
strictly vanish. Thus there will also be some constraint-violating processes. 
Second, while lattice models are often good approximations for static properties, 
for dynamics of systems in which the geometry and/or bonding is believed to 
play an important role, it is dangerous to appeal to universality arguments for 
the validity of lattice models. This is especially true if the length scales of 
the important processes never get very large, as is believed to be the case for 
glasses [3]. 

Both of these issues merit much further attention. In particular, one should 
at a minimum consider the effects of a low rate of constraint-violating particle 
motions within the lattice model framework. In practice, the neglect of certain 
processes may be a good approximation on a wide range of time scales. For 
example, if there are chemical bonds with energies much larger than fc^T, then 
the ratio of the rates of processes between, e.g., those that require breaking of 
two bonds and those that require breaking three or more, may be very small. A 
— and arguably the — fundamental problem of dynamics near glass transitions 
is the large and rapidly-onsetting increases in the effective activation barriers. 
Thus looking for an explanation that involves a rapid crossover from domination 
by a relatively fast process to freezing out of this and the ensuing domination 
of relaxation by a much slower process, is surely a productive direction. 

The issue of continuum versus lattice models is a trickier one. If geometrical 
constraints really dominate near glass transitions, then one needs to consider 
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how to model the interplay between energetic and entropic barriers, even if the 
former only arise from repulsive interactions that are not entirely hard-core. 
In principle, from simulations or from computations of more realistic contin- 
uum models of molecular interactions, one may be able to extract effective rates 
for various local processes (see [12] for an investigation in this direction). If 
temperature changes dominate over density changes [68] the effective density 
parameter in a lattice approximation will not be simply related to the actual 
density because of soft-core effects. Nevertheless, if free-volume ideas of glass 
transitions are relevant, one should be able to work in terms of some effective 
density — with "vacancies" , or at least mobile ones, representing the free vol- 
ume. Because of local expansion and contraction, however, the free volume is 
in general not conserved. One might expect the particle conservation to play a 
crucial role in KA models. In fact, this is not the case, as the non-conserving 
Frcdrickson- Anderson model shows very similar behavior to KA models. Pre- 
sumably, a better approximation would conserve (effective) density most of the 
time, but not always. 

We leave these issues and any attempts at direct comparisons with experi- 
ments, for future papers. 
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